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FOREWORD 


This  study  was  conducted  by  the  Department  of  Oceanography  at  Texas  A&M 
University,  Reference  85-2-T.  The  research  was  conducted  through  Texas  A&M 
Research  Foundation  Project  4667. 

This  study  was  motivated  by  the  fact  that  existing  storm  surge  models 
which  are  generally  restricted  to  a  limited  reach  of  the  continental  shelf 
do  not  adequately  simulate  the  forerunner  surge.  In  the  Gulf  of  Mexico  the 
forerunner  stage  of  a  surge  event  (which  can  have  an  amplitude  as  large  as 
1  m)  can  precede  the  peak  surge  by  more  than  24  hours. 

It  is  common  practice,  in  employing  coastal  surge  models  for  design 
purposes,  to  start  computations  with  an  initial  non-zero  constant  water 
level  field  over  the  shelf  in  recognition  of  the  influence  of  forerunners. 
This  procedure  makes  no  allowance  for  initial  flows  or  gradient  thereof 
which  must  exist  during  the  forerunner  stage.  On  the  other  hand,  starting 
local  coastal  models  well  in  advance  of  the  arrival  of  the  hurricane  at  the 
shelf  break  cannot  simulate  the  forerunner  surge  since  the  latter  is  ex¬ 
cited  by  mechanisms  acting  over  the  adjacent  deep  basin  seaward  of  the 
coastal  rone  being  modeled. 

The  objectives  of  this  research  were  to  determine  for  the  Gulf  of 
Mexico  the  space  and  time  scales  of  forerunners,  examine  possible  excita¬ 
tion  processes,  and  evaluate  the  role  of  baroclinic  modes  in  surge  events 
at  the  c<>ast  and  shelf  break.  A  numerical  approach  was  used  in  this  study 
which  modeled  hurricane-forced  responses  of  the  barotropic  and  first  baro¬ 
clinic  modes  over  a  computing  domain  representing  the  Northwest  Cayman  Sea 
and  Gulf  o*  Mexico  The  inclusion  of  the  baroclinic  response  adds  a  degree 
of  realism  acking  in  previous  studies  on  hurricane  effects  in  the  Gulf  of 
Mexico. 
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It  should  be  stressed  that  the  model  includes  only  wind  and  atmospheric 
pressure  forcing.  Tidal  forcing  was  not  included.  Moreover,  the  model  is 
quasi-linear  (only  the  bottom  and  interface  stresses  are  nonlinear) .  The 
model  w.^s  verified  by  simulating  the  Gulf's  response  to  hurricanes  Carla  and 
Allen.  Although  the  focus  of  the  study  was  on  forerunner  surge  evei._-,  the 
results  also  provide  insight  into  fundamental  concepts  concerning  the  large- 
scale,  low-frequency  free  modes  in  the  Gulf  of  Mexico  and  Cayman  Sea. 

This  study  also  constitutes  the  doctoral  dissertation  research  of  one 
of  the  authors,  Mahunnop  Bunpapong. 
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CHAPTER  I 


INTRODUCTION 


1.  Background 

Redfield  and  Miller  (1957)  considered  the  changes  in  water  level 
associated  with  hurricanes  as  consisting  of  three  successive  stages; 
the  forerunner  surge,  the  hurricane  surge  and  the  resurgence.  The 
forerunner  surge  is  the  gradual  rise  in  water  level  along  the  coast 
which  precedes  the  arrival  of  the  hurricane.  It  occurs  while  the 
storm  center  is  at  a  great  distance  from  the  coast,  irrespective  of 
whether  the  storm  reaches  the  point  of  observation.  The  gradual 
buildup  of  water  level  may  reach  one  mete”.  This  may  seem 
unimportant  in  a  destructive  sense  when  compared  to  the  hurricane 
surge  itself,  but  it  is  an  important  initial  boundary  condition  for 
local  storm  surge  models. 

The  salient  feature  of  forerunners  is  that  local  atmospheric 
forcing  is  not  required  tor  this  phase  of  the  surge.  Cline  (1920) 
reported  the  existence  of  forerunner  surges  in  the  Gulf  of  Mexico. 

He  analyzed  water  levels  and  winds  associated  with  hurricanes  from 
the  1900-1919  period  and  found  that  some  storms  clearly  produced 
increasing  sea  level  within  24  h  of  the  disturbance  entering  the 
Gulf.  Cline  postulated  that  forerunners  are  related  to  swell 
generated  by  Gulf  hurricanes.  Based  on  a  simple  wind  model,  Cline 
showed  that  the  storm  generated  swell  (forerunner)  would  be 
significant  only  in  the  direction  of  the  hurricane  motion.  The  data 
presented  by  Cline,  however,  suggest  the  nearly  simultaneous 
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appearance  of  forerunners  over  a  long  stretch  of  the  United  States 
Gulf  coast. 

Forerunners  in  a  basin  like  the  Gulf  of  Mexico  may  be  explained 
in  terms  of  the  large  scale  barotropic  normal  gravity  modes  of  the 
Gulf.  Even  though  the  Gulf  is  relatively  small,  the  size  of 
hurricanes  still  gives  a  mismatch  of  forcing  and  response  scales. 
However,  normal  mode  excitation  is  not  limited  tc  open  ocean 
atmospheric  forcing  in  the  Gulf  of  Mexico.  Platzman  (1972)  obtained 
a  long-period  Helmholtz  mode  of  21.2  h  in  a  numerical  study  of  the 
Gulf.  Reid  and  Whitaker  (1981)  also  obtained  a  Helmholtz  mode  for 
the  Gulf,  but  with  a  period  of  28.5  h.  This  mode  is  characterized  by 
nearly  uniform  phase  and  amplitude,  except  near  the  ports. *  This  is 
indicative  of  a  co-osci Hating  port-driven  Gulf -Caribbean  system. 
Conceivably,  some  hurricanes  could  provide  port  forcing  to  olicit 
this  long  period  motion. 

Another  possible  large  scale  barotropic  response  of  the 
hurricane  forcing  typically  has  relatively  small  surface  expressions 
but  large  vorticity  (circulation).  These  are  the  vorticity  modes 
generated  by  planetary  or  topographic  vortex  stretching  and  they  have 
been  observed  in  enclosed  basins.  Saylor  et.  al  (1980)  found  that 
the  observed  oscillation  in  Southern  Lake  Michigan  with  a  period  of  i 
days  was  a  lake-wide  barotropic  topographic  vortex  mode. 

Calculations  show  that  if  the  response  of  the  Gulf  of  Mexico  to  a 
hurricane  was  a  Gulf-wide  barotropic  vortex  mode  it  would  have  a 
period  of  about  5  days.  Kullar  and  Prankignoul  (1981)  however, 
showed  that  realistic  topography  causes  the  vorticity  trades  to  have 

*  The  tern  “porta"  denotes  the  model  representation  of  the  Yucatan  and 
Florida  Straits. 


smaller  space  scales  compared  to  the  gravity  modes. 

The  development  of  forerunners  in  the  Gulf  might  be  caused  by 
the  strong  baroclinic  response  of  the  sea  to  hurricane  winds. 
Stevenson  and  Armstrong  (1965)  and  Leipper  (1967)  presented 
descriptions  of  hurricane  induced  upwelling  in  the  Gulf  and  provided 
qualitative  depictions  of  the  causative  mechanisms.  Oceanic 
baroclinic  responses  to  hurricane  forcing  have  been  the  subject  of 
numerous  theoretical  ar.C  numerical  investigations  (Ichiyc,  1955).  In 
addition,  Veronis  and  Stonmel  (1956)  speculated  that  a  two- layer 
ocean  with  constant  depth  would  exhibit  a  baroclinic  response  to  a 
distant  storm.  However,  changes  in  the  interface  elevation  would  be 
small  provided  that  there  were  no  other  effects  such  as  resonance. 

The  only  noticeable  effect  would  be  changes  in  the  free  surface 
anomaly  contributed  by  barotropic  Rossby  waves. 

The  single  study  on  forerunner  surges  in  the  Gulf  by  Cline 
(1920)  leaves  several  fundamental  questions  unanswered.  It  is  not 
clear  that  all  hurricanes  traversing  the  Gulf  generate  forerunners. 

If  not,  under  which  conditions  are  forerunners  excited?  What  are  the 
time  and  space  scales  of  forerunners?  Conflicting  evidence  suggests 
both  local  avid  Gulf  wide  occurrences.  Finally,  are  fc^er'nners 
barotropic  or  baroclinic  in  origin?  To  date,  only  the  barotropic 
response  on  the  perimeter  of  the  Gulf  to  hurricane  forcing  has  bean 
sufficiently  studied. 

As  already  pointed  out,  the  forerunner  questions  are  important 
in  establishing  initial  conditions  in  hurricana  surge  prediction  at 
the  coastline.  Most  hurricane  surge  models  employ  limited  area 


domains,  namely  a  section  of  the  continental  shelf  extending  from 
shore  to  the  shelf  break  at  about  the  200  to  depth  contour  and 
extending  several  hundred  kilometers  along  shore  on  either  side  of 
the  coastal  location  for  which  surge  prediction  is  sought.  A  second 
problem  with  such  models  is  that  of  specifying  appropriate  boundary 
conditions  at  the  open  boundaries,  particularly  along  the  seaward 
shelf  break  positions.  A  common  seaward  boundary  condition  is  to  set 
the  water  level  along  the  shelf  break  equal  to  the  static 
barometrically-induced  anomaly  appropriate  to  the  position  relative 
to  the  hurricane  center.  However,  this  ignores  any  dynamically- 
induced  water  level  anomaly  by  the  storm  in  the  deep  region  of  the 
sea  or  Gulf. 

2.  Objectives 

There  are  three  objectives  of  this  research  effort.  First,  the 
investigation  is  directed  toward  establishing  the  cause  as  well  as 
the  time  and  space  scales  cf  forerunners  in  the  Gulf  cf  Mexico. 
Second,  the  relation  between  hurricane  generated  barotropic  and 
baroclinic  modes  and  forerunners  is  sought,  nurd,  the  clarification 
of  conditions  at  the  shelf  break  is  addressed.  All  of  these 
objectives  will  be  carried  out  within  the  framework  of  a  basin  wide 
model  which  includes  the  Gulf  of  Mexico  and  the  Cayman  Sea. 


3.  Procedure 


As  a  result  of  uncertainty  in  the  type  of  response,  it  is 
prudent  to  include  the  effect  of  density  stratification  in  the 
formulation  of  the  problem.  To  allow  for  the  first  baroclinic  mode 
in  the  solutions,  a  two- layer  model  with  variable  depth  was  chosen. 
An  analytic  approach  to  ths  problem  of  hurricane  forcing  in  a  basin 
with  realistic  bathymetry  and  shape  is  not  possible.  Therefore,  a 
numerical  approach  was  used  in  this  study.  An  existing  linaar 
numerical  model  for  astronomical  tides  in  the  Gulf  of  Mexico  ,GOMT, 
(c.f.  Reid  and  Whitaker,  1981)  was  adapted  to  include  both  the 
barotropic  and  baroclinic  computations  and  atmospheric  forcing. 
Information  obtained  from  the  numerical  model,  which  is  usually 
unavailable  from  the  recorded  data,  allows  a  detailed  investigation 
of  the  dynamics  and  other  characteristics  of  the  response.  Another 
important  advantage  of  the  numerical  model  is  its  predictive 
capability. 

Instead  of  using  the  primitive  equations  describing  the  two- 
layer  system,  normal  mode  equations  were  derived  using  a 
generalization  of  the  method  employed  by  Veronis  and  Stomnel  <1956). 
There  are  several  advantages  in  working  with  the  modal  form  of  the 
equations.  In  this  system  the  dynamics  and  energetics  of  each  mode 
can  be  examined  separately.  Interaction  between  modes,  due  to 
coupling  caused  by  varying  depth,  can  be  investigated  in  terms  of 
energy  transfer.  Open  boundary  conditions  are  facilitated  in  terms 
of  the  modes,  particularly  if  outward  radiation  of  wave  energy  is 
allowed.  In  addition,  the  layer  variables  such  as  the  free  surface 
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CHAPTER  II 


THEORETICAL  AND  NUMERICAL  MODEL 

1.  Theoretical  model 


a)  Basic  equations 


The  vertically  integrated  momentum  and  mass  conservation 
equations  for  quasi-hydrostatic  (large  scale)  disturbances  in  a  two- 
layer  variable  depth  basin  are,  for  the  upper  layer. 


aMj/at  +  fkxM1  +  p1gH1V(h1+h2)  =  P1# 


Pi  ahj/at  ♦  =  o. 


and  for  the  lower  layer, 


*>  -» 


3M2/3t  +  fkxM2  "*■  gH2^(p^h^,*p2h2 )  ®  F2, 


p2  ah?/at  ♦  v-m2  =  o, 


(1) 

(2) 

(3) 

(4) 


where  M  is  the  mass  transport  per  unit  width,  f  is  the  Coriolis 
parameter,  k  is  the  vertical  unit  vector,  p  is  the  water  density, 
g  is  the  gravitational  acceleration,  H  is  the  mean  depth,  H*h  is  the 
instantaneous  depth,  F  is  the  external  forcing  and  dissipation 
defined  as: 


P1  *  Ts  “  Ti  ~HlV5?a' 

f2bVV  h2*V 


(5) 
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Here  r  is  the  stress  vector  where  subscripts  s,  i,  and  b  stand  for 
surface,  interface,  and  bottom,  respectively.  The  atmospheric 
pressure  at  the  sea  surface  is  Pa. 

The  normal  mode  form  of  the  equations  can  be  derived  from  these 
primitive  equations  by  a  method  similar  to  that  employed  by  Veronis 
and  Stommel  (1956).  To  transform  (1)  and  (2)  into  normal  mode  form, 
we  multiply  (1)  and  (2)  by  a  and  (3)  and  (4)  by  (5  and  add  the 
corresponding  momentum  and  mass  conservation  equations,  respectively, 
to  obtain 


3M/3t  ♦  fkxM  ♦  g{V[aH1(p1h1+p1h2)>0H2(p1h1+p2h2)] 

“ [  ( P^h^>p^h2 ) V(aH^)  +  ( p^h^^p2h2) V( ()H2)  ]}  =  G, 

d*/3t  ♦  ♦  0V-M2  =  0, 


where 


H  m  aM^  ♦  0H2, 

*  ts  ap^  ♦  0P2h2, 

■+  -*  -* 

G  s  aF^  ♦  (J?2. 


(6) 

(7) 


(8) 


Note  that  a  and  0  are  non-dimensional  and,  for  the  case  of  variable 
depth,  may  depend  upon  x  and  y  (this  is  the  generalization  of  the 
Veronis  and  Stommel  analysis). 

The  constraint  imposed  on  a  and  0  to  make  the  elevation 
anomalies  in  (6)  and  (7)  proportional,  is 


cH1(pihi»p1h2)  ♦  0H2(p1h1*p2h2)  ■  r#. 


(9) 
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where  P,  a  factor  of  proportionality,  is  an  equivalent  depth  to  be 
determined.  Since  $  =  ap^h^  ♦  fiP2tl2>  then  O)  will  be  valid  for  all 
combinations  of  h^  and  h2  if 


aH^  *  0H2  a 

rc,  and 

(10a) 

aHl ( P\/ ^2 )  +  ^2  = 

P0. 

(10b) 

It  can  be  readily  shown  that  (6)  and 

-*  -*  -*  -» 

(7)  can  be  written 

in  the  form 

3M/3t  ♦  fkxM  ♦  grVs6  -  B  = 

G, 

(11) 

« 

3*/3t  ♦  V-M  -  C  = 

0  , 

(12) 

where 


B  e  gHpjhjVt  ♦  P2h2?0)  ,  and 

•*  •* 

C  s  M^Va  ♦  M2*V0. 


(13) 


For  the  case  of  constant  layer  depth  H2#  the  a  and  0  are 
*• 

constant  and  B,  C  vanish.  In  general  these  terms  produce  coupling  of 
the  inodes  in  the  presence  of  bottom  topography. 

Eqs .  (10a,b)  can  be  written  in  matrix  form  as 


'  Hrr 

“2* 

m  ■“ 

a 

V 

_(p1/p2)H1 

H2-r 

.0. 

.0. 

(14) 


and  hence  the  eigenvalues  of  r  are  the  roots  of  the  characteristic 
equation 


9 


r*  -  or  +  «h1h2 


*  o  , 


(15 


where  D  is  the  total  depth  at  rest  (H^+H2) ,  and  <  is  the  relative 
density  difference,  (p 2“pl^p2* 

The  two  roots  of  (15),  which  correspond  to  the  equivalent  depths 
for  external  and  Internal  modes,  are 


re  =  D[1-«(H1B2/D2)  ♦  0(«2)]  ,  and 

rL  =  «(h1h2/d)[i*.(h1h2/d2)  ♦  o(«2)]. 


(16: 


Since  H1H2/D2 


<  1  and  t  «  1  for  the  general  ocean  basin, 


then 


re  =  D  ,  and 


rL  a  «(h1h2/d) 

are  the  equivalent  depths  for  external  and  internal  modes  as  obtained 
by  Veronis  and  Stommel  (1956).  In  general  these  may  depend  on  x,y. 
The  ratio  of  a  and  (3  can  be  determined  from  (10a,b)  as 

a/a  =  r  -  Hj/Hj  »  (pj/pjMj/r  -  h2.  (is: 

Substituting  (17)  in  (13),  the  ratio  (a/St)  for  each  mode  is 
(a/a)e  -  1  -  fUij/D)  -  «2(h|h2/D3)  ,  and 

(19) 

(a/nJi  -  -H1/H2  ♦  tUij/D)  ♦  «2(HjH2/D3). 


Note  that  ae  and  j?e  are  of  like  sign,  but  and  0^  are  of  opposite 
sign.  ce,  ac,  a±  are  chosen  to  be  positive  and  ai  to  bo  negative. 


Veronis  and  Stommel  (1956)  took  ae  *  *  1  for  a  constant  depth 

basin.  However,  the  presence  of  coupling  terms  B  and  C  in  (11)  and 
(12)  due  to  varying  depth  in  the  present  study  precludes  an  arbitrary 
choice  of  c.  Notice  that  both  a  and  0  are  functions  of  x  and  y. 

The  procedure  used  to  determine  the  individual  values  of  a  and  0 
for  both  modes  is  based  on  energy  considerations.  In  essence  the 
energy  equations  derived  from  the  primitive  equations  and  the  normal 
mode  equations  must  be  consistent. 

The  two  energy  equations  formed  from  the  primitive  equations  and 
the  normal  mode  equations,  respectively,  are 

3(Ek+Ep)/at  ♦  V-J  =  S,  (20 

where 

Ek  a  *  2p^h^h2  ♦  (>2^2")  > 

Ep  =  JKM^/P^)  ♦  (M2/p2H2), 

(21 

J  c  g[(h1>h2)Mi  ♦  1/ »2 ^ Pihj^p2h2 )M2 ] , 


S  ~  H^*F^/p^H^  ♦  M2-F2/p2H2. 

and 

d(Ek>Ep)/at  ♦  V-J  -  T  ■  S  .  (22 

Equation  (22)  contains  not  only  the  additional  term  T,  but  the 

expressions  for  the  kinetic  energy  per  unit  area,  E^,  the  potential 

*♦ 

energy  per  unit  area,  Ep,  the  energy  flux  per  unit  width,  J,  and  the 
net  energy  supply  per  unit  width  per  unit  time,  S,  are  also  different 


from  those  given  in  (21).  The  added  term,  T,  which  defines  the 


energy  transfer  between  modes,  and  the  other  terms  in  (22)  are 
defined  as  follows: 

T  =  l/p2  [  <  V  Vre  >  ♦  « <  Mi '  Vri )  +g  <  *eV  *  )  ] .  (23) 


Ek  =  l/p2[(M|/re)  ♦  «(Mf/ri)],  (24a) 

Ep  =  l/p2gU|  ♦  «42}  •  (24b) 

J  =  g/p2(^eM0  ♦  i$^M^),  (24c) 

S  =  l/p2[(Se*Ge/reH«(Mi.Gi/ri)],  (24d) 

where  Ge  and  sure  defined  in  (31). 


If  (20)  and  (22)  are  to  be  consistent,  their  corresponding  terms 
must  be  equal  and  T  must  be  zero.  Inserting  (8a)  to  (8c)  in  (24a)  to 
(24d)  it  can  be  shown  that  the  necessary  conditions  that  Ek  as  well 

•*  *4 

as  S  in  (20)  and  (22)  are  consistent  for  arbitrary  hj_,  h2,  M1#  and  M2 


are  that 

ae2/r@  ♦  «ai2/ri  =  l/H1(p2/p1),  (25a) 

0e2/re  ♦  *^i2/ri  58  1/K2*  (25b) 

oe5e/re  ♦  .o^i/ri  =  0.  (25c) 

It  is  also  necessary,  for  consistency  in  Ep  and  J,  that 

oe2  ♦  •oi2  ■  p1/p2,  (26a) 

fle2  ♦  «i5i2  -  1  »  (26b) 


oJL  ♦ 


»  1  . 


(26c) 


The  requirement  that  T  must  be  zero  can  be  examined  by  substituting 
(13)  in  (23)  and  rearranging  to  yield, 

P2h2*V7(0e2*«%2)-  <27 ) 


Using  conditions  (26a)  to  (26c),  it  can  be  shown  that  T  vanishes 
provided  that  the  ratio  p2/ p^  is  constant.  Therefore,  (25)  and  (26) 
assure  that  (20)  and  (22)  are  entirely  consistent  with  T  equal  to 
zero  and  all  other  terms  are  identical.  These  two  sets  of  equations 
will  also  be  used  to  determine  a  and  0. 


Prom  (25a), (25b)  and  (26a), (26b)  we  find 


and 


Ge2  *  P2re/^>lHl^Hl"ri)/(re'ri)^ 

0e2  *  VH2[(H2-ri)/(re-ri)], 


(28) 


°i2  ■  p2ri/«P1H1[(re-H1)/(re-ri)], 
Oi2  »  ri/«H2[(re-H2)/(re-ri)]. 


(29) 


It  is  not  difficult  to  show  that  the  ratios  (0/o)e  and  (0/a)^ 
obtained  from  (28)  and  (29)  are  the  same  as  those  in  (19).  If  we 
substitute  (16)  in  (28)  and  (29),  we  obtain  the  relation  for  a  and 
0, accurate  to  order  «,  in  terms  of  H^,  H2  and  D  as  follows: 

ce  *  1  ♦  i*  ~  J«(H2/D)2  ♦  0 («2),  (30a) 

0e  -  1-  J.(HX/D)2  ♦  0(«2), 


(30b) 


aL  =  H2/D{1^«  ♦  «[(H1H2/D2)-J(H1/D)]  ♦  0(«2)}  (30c) 
fiL  =»  -H1/D{l+f[(H1H2/D2)-J(H2/D)]  ♦  0(«2)}  (30d) 


Note  that  positive  roots  are  chose-i  except  for  Using  a  and  0 
from  (30)  it  can  be  verified  that  (25c)  and  (26c)  are  correct  to  the 
order  «. 


b)  Normal  Mode  Equations 


To  express  Eqs.(ll)  and  (12)  with  all  dependent  variables  in 
terms  of  modes,  (8)  is  rewritten  for  each  mode  as: 

Me  =  ♦  0eM2 


*e  3  aeplhl  *  ^ep2h2 


♦i  *  °iplhl  *  ^ip2h2 


Ge  “  aeFl  *  0ep2 


Gi  3  <*iFl  ♦  0ip2' 


It  readily  can  be  shown  from  (31)  that 


plhl  *  r(Vi  "  5iV 


p2h2  -  ♦  <4#e) 


M1  3  r^eKi  ’  W 


M2  s  r ( ♦  a^Me) , 

where  r  =  (a^g  -  ae0i)-i.  Using  (30a)  to  (30d),  the  factor  r  is 
(P]/P2)1/,Z  to  or<ier  «  and 


h2  =  [$  ♦  (H1/D)^e](p1p2) 

h2  =  [<«H2/D  -  1)*  ♦  (H2/D)^e](p1p2) 


Mx  =  (H1/D)Me  +  Mi, 


4' 


(33) 


4  4 


M2  «  (H2/D)Me  -  Mi. 


Using  (31)  and  (32)  the  coupling  terms  defined  by  (13)  can  be 
approximated  by. 


De  a  r9ree*i7(Hl/D)' 

Bi  a  -gri^Hj/D), 

Ce  a  r.Mi-VCHj/D), 
q  a  -Mg-VtHi/D). 


(34) 


These  forms  show  clearly  that  the  external  mode  is  influenced  by 
the  internal  mode  and  vice  versa  when  H2/D  is  variable. 

It  can  be  shown  by  substituting  (34)  in  (23)  that  the  energy 
transfer  between  modes,  T,  is  aero  as  required. 

Finally,  the  forcing  terms  G  can  be  approximate!  by  using  (S), 
(8),  and  (31)  as: 


IS 


Ts  “  Tb  “  D^pa  4  °(*)* 


(35] 


Gl  =  (H2/D>fs  ♦  (H1/D)rt>  -  TL  *  0 («). 


counterparts  of  M, 

0,  r  and 

pa 

are  defined 

-♦ 

-* 

-> 

Qe 

■ 

Me/P. 

Qi 

■ 

M^p, 

^e 

S 

Vp' 

*i 

B 

^i/p* 

-* 

-♦ 

Ts 

Id 

rs/p, 

Ti 

■ 

r  i/p. 

-♦ 

-• 

Tb 

5 

t^/p, 

-vpa/p 

■ 

gVb, 

(36) 


The  variable  b  in  (36)  is  the  barometric  pressure  deficit  expressed 
in  terms  of  an  equivalent  head  of  water  (sometimes  referred  to  as  the 
inverted  barometric  effect). 

The  normal  mode  equations  (11)  and  (12),  therefore,  take  the 

form 


3Qe/3t  ♦  fkxQ@  ♦  gDV*e  -  «gJtyiV(H1/D)  =  (Ts  -  Tb  ♦  gDVb), 


(37: 


3ve/3t  ♦  V-Qe  -  *  0  , 


and 


aQj/dt  ♦  fkxQi  ♦  griV^i  ♦  gri^eV(H1/D)  »  (H2/D)Ts  ♦  (Hj/D)Tb  -  xL 


i » 


(3s: 


3Vj/3t  ♦  r  Qi  ♦  Qe*V(H1/D>  »  0 


where  r,  is  given  in  (17). 


2.  numerical  model 


a)  Grid  system 

A  space-staggered  computational  mesh  is  employed  in  the 
numerical  analogs  of  the  normal  mode  equations  (37)  and  (38).  The 
grid  spacing  is  taken  as  15’  x  15*  in  latitude  (A$)  and  longitude 
(AX) .  The  modal  transport  per  unit  width,  represented  by  U,  V  for 
east  and  north  components,  respectively,  is  located  at  the  mid  points 
of  the  appropriate  sides  of  each  grid  block.  The  modal  water  level 
anomaly,  4>,  is  defined  at  the  middle  of  each  grid  block.  Figure  1 
illustrates  locations  of  these  variables  whose  positions  are 
identified  by  I,J  and  time  is  indexed  by  n.  The  grid  spacing  is 
taken  as  the  distance  between  the  same  variables,  thus  consecutive  U, 
4/  values  are  at  half  increments  of  I  and  consecutive  v,  ^  values  are 
at  half  increments  of  J.  The  depth  is  stored  at  every  variate 
location. 

The  grid  system  used  in  this  study  is  the  same  as  that  used  in 
the  Gulf  of  Mexico  Tidal  Model  (Reid  and  Whitaker,  1981).  It 
nominally  covers  70  x  50  grid  blocks  for  the  Gulf  of  Mexico  and  the 
northwest  portion  of  the  Caribbean  Sea  (Cayman  Sea)  as  shown  in  Fig. 
2. 

The  depths  of  the  Gulf  of  Mexico  were  digitized  from  bathymetry 
charts  prepared  by  Buffler  et  al.  (1984).  The  Cayman  Sea  depths  were 
taken  from  available  hydrographic  charts.  Depth  over  most  of  the 
Cayman,  especially  near  the  open  boundaries  to  the  east  and  south, 
was  taken  as  4000  a.  Figure  3  shows  the  computer  plotted  contours  of 
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vv — H 


Computing  lattice  showing  locations  where  U,  V,  ^  are  evaluated.  Positive  U  is  toward  the 
increasing  index  I).  Positive  V  is  to  the  north  (increasing  index  J). 


system  for  the  Gulf  of  Mexico  and  the  Cayman  Sea.  The  grid  Increments 


gjfjfjggjgjgj 


coastline  and  sample  contours  of  the  digitized  bathymetry.  Depths  in  meters 


the  digitized  depth  field.  Special  care  was  observed  when  comparing 
the  digitized  depth  fields  to  the  bath;T>ef.ry  chart,  especially  in  the 
shelf  break  regions  and  modifications  were  made  where  necessary. 

b)  Numerical  integration  scheme 


The  multioperational  alternating  direction  implicit  algorithm 
developed  by  Leendertse  (1967)  was  adopted  for  time  integration  of 
the  finite  difference  equations.  The  following  notation  is  used  in 
the  discussion.  The  spatial  average  of  a  field  variable  X  is 
written  as 

Xn(I,J)  *  ■i[X(I-^#J-i)+X(I-^,J-^)+X(I-J,J^)+X(J-^,j4)}»  (39) 

where  Xn(I,J)  s  X(X0+lAX,4>0+JM>,t0+nAt) .  Time  and  space  derivatives 
are  depicted  by  the  standard  centered  differences, 


ax/at  -  (i/At)[xnn(i,j)-xn_1(i,j)3, 
ax/ax  =  (cos$/&x)[xnu4rJ)-xn(i-^,j)]. 


(40 


k  spherical  coordinate  system  is  employed  in  representing  the 
gradient  terms  on  a  level  surface,  i.e., 


VX  =  (l/a0( J) )  {aX/ 34*  a  ♦  9(J)3X/3X  b} 


(41 


where  X  is  any  scalar  field  variable,  a  is  a  unit  vector  along  lines 

*4 

of  constant  <J>,  b  is  a  unit  vector  along  lines  of  constant  X,  and 
0(J)  »  cos(4>o+JM),  where  is  the  reference  latitude  (18*N). 
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The  cycle  of  calculation  is  separated  into  two  operations. 

During  the  first-half  cycle,  at  odd  time  steps,  and  U  are  computed 
implicitly  along  lines  of  constant  latitude,  followed  by  an  explicit 
computation  of  the  V  field.  For  the  second  half -cycle,  the 
computations  proceed  along  lines  of  constant  longitude  with  p,  V 
updated  implicitly  and  U  computed  explicitly.  In  each  half -cycle, 
the  external  mode  computations  are  executed  first  and,  after 
completing  the  entire  computing  domain,  are  repeated  for  the  internal 
mode. 

The  implicit  formulation  of  the  finite  difference  analogs  of  the 
external  mode  momentum  and  mass  conservation  equations,  respectively, 
are.  for  odd  time  steps, 

,  n+1  ,  n+1  i  n+1  . 

~VI4''me(I~?'J)  4  ue(I'J>  *  ?x(l4Z'J)*e(l42'J) 

*  Ug( I, J)+2Atf (J)Vg(I, J)  ♦  AtFfi*  XL,  and  (42) 

n+1  i  n+1  n+1  ■. 

-i.x(J)Ue(l4,J>  +  *e(I,J>+  ux(J)Ue(I+£,J> 

=  *J(I,J)  -  uyfiyv£(I,J)  +  {i#  (43) 

and  for  even  time  steps 

-7y(i,j4)4u.j4>9(J4>  4  Vgti, j)  +  7ydfj4)v^J4)d(J4i) 

*  v£(I,J)  -  Atf(J)u£(I,J)  +  Xi+  AtFe,  (44) 

,  n+1  ,  n+1  ,  n+1  , 

-Wy(j4)Ve(l,j-^)  +  *e(I,J)  ♦  uy(j4)Ve(I,j4) 

“  '  vx5xue(I*J)  4  *i*  (45) 
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where 


7X<I,J) 

=  (At/AX) {gD(I,J)/afi(J)}, 

7y(I»J) 

=  (At/A*){gD(I,j)/a0(J)}, 

UX(J) 

=  (At/AX) {l/a«(J)}, 

»y(J) 

=  (At/A*){l/afl(J)},- 

SKuJ(I,J) 

=  UJ(I+J,J)  -  u{j(I-$,J), 

8yVj(l,j) 

=  VJ(I#J+J)«(J+J>  -  vJ(I,J-J)fi(J-J), 

a  =  radius  of  the  Earth, 
f(J)  =  2fisin(4>0+JA$) , 

F  *  forcing  and  friction  terms, 

Xi$  3  coupling  terms. 

The  forcing  and  coupling  terms  will  be  discussed  in  the  next  section. 
The  coupling  terms  in  (42)  and  (43)  or  (44)  and  (45)  are  of  opposite 
mode  relative  to  the  other  terms. 

Upon  replacing  the  total  depth,  D,  by  the  equivalent  depth, 
eHiH2/D,  interchanging  modes  of  all  field  variables  and  coupling 
terms,  and  using  the  proper  forcing  terms,  (42)  through  (45)  are  also 
representations  of  the  implicit  formulation  for  the  internal  mode 
computation. 

The  explicit  coding  of  the  external  mode  momentum  equation  at 
odd  time  steps  is 
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(46) 


VI, J)  =  VJ(I,J)  -  jAt  f(J){Ug(I,J)  +  ug(I,J)} 
-7y(I/J)&y^(I,J)  +  Xi  +  AtFg. 

The  even  time  step  counterpart  of  (46)  is 

Ug(i,j)  =  uJ(i,j)  -  JAt  f(j){yi,j)  +  vg(i,j)} 

-7x(I,J)6x^(I,J)  +  Xi  +  At  Fg.  (47) 

Likewise,  the  explicit  coding,  at  odd  and  even  time  steps, 
respectively,  for  the  internal  mode  computation  is  obtained  by 
replacing  the  total  depth  by  the  equivalent  depth,  interchanging 
modes  of  all  field  variables  and  coupling  terms,  and  employing  the 
proper  forcing  terms  in  (46)  and  (47). 

Eqs.  (42)  and  (43)  or  (44)  and  (45)  form  a  system  of  linear 
algebraic  equations  in  the  collective  I  or  J,  depending  upon  time 
step,  of  4/  and  either  U  or  V  at  time  level  n+1.  The  coefficient 
matrix  is  tridiagonal  for  which  there  exists  a  double  sweep  solution 
algorithm  for  inversion,  provided  that  boundary  conditions  on  U  or  V 
or  some  combination  of  conditions  on  U  or  V  and  \p  are  supplied  at 
each  end  of  the  array  of  variables. 

c)  Surface,  interface,  and  bottom  stress 

The  forcing  term  F  in  (42)  and  (44)  consists  of  the  surface 
stress,  the  bottom  stress  and  either  the  atmospheric  pressure  force 
due  to  a  surface  pressure  deficit  for  the  external  mode  forcing  or 
the  interface  stress  for  the  internal  mode  forcing,  respectively. 

That  is 
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Fe  =  Ts  -  Tb  +  gDVb,  or 

(48) 

H,  -♦  H,  ->  -» 

Fi  =  <Ef  >Ts  +  <13  >Tb  -  Ti* 

The  stress  terms  are  presented  in  the  form 


T  =  *|w|w.  (49) 

*4 

For  the  surface  stress,  W  is  the  wind  speed  at  an  elevation  of 
10  m  above  the  water  surface.  Reid  and  Bodine  (1968)  considered  *c  as 
a  function  of  wind  speed  in  the  form 


K  =  K}. 


for  |w|  S  Wc 


(50) 

k  =  kx  +  k2(1-Wc/|w|)2  for  | W |  ^  Wc 

where  k1  and  <2  are  taken  as  1.1  x  10-6  and  2.5  x  10-6,  respectively, 
and  Wc  is  a  critical  wind  speed  which  is  taken  as  7.0  m/s.  The 
coefficient  k  is  related  to  the  drag  coefficient,  CD  by  the  relation 


*  =  (^a^  ^w^D 


(51) 


where  pa  is  air  density  and  pw  is  water  density. 

For  large  wind  speed,  k  approaches  the  limiting  value  of  3.6  x 
10~6  which  corresponds  to  a  drag  coefficient  of  about  3.0  x  10~3  if 
the  density  ratio  between  air  and  water  is  assumed  to  be  1.2  x  IQ-3. 
Equation  (50)  was  used  by  Wanstrath  (1975)  in  his  simulation  of  storm 
surge  in  transformed  coordinates  while  Miyasaki  (1963)  used  a 
constant  3.2  x  10~6  for  k  in  his  computation  of  storm  surge  for 
hurricane  Carla.  The  choice  of  k  for  intense  winds  (|W|£  50  m/s)  is 


25 


controversial.  However,  taking  into  account  that  wind  speeds 
associated  with  hurricanes  are  not  steady  or  uniform,  *  is  taken  as  a 
function  of  wind  speed. 

The  surface  wind  stresses  are  computed  at  every  time  step  using 
linear  interpolated  positions  of  the  hurricane  center  which  are  given 
at  6  h  intervals,  except  for  the  case  of  hurricane  Carla  simulation 
which  will  be  discussed  later.  A  constant  inflow  angle  of  20°  is 
used  to  rotate  the  surface  stress  vector  before  the  components  are 
computed . 

■4 

For  the  bottom  stress,  W  is  the  depth  averaged  lower  layer 
current  velocity,  i.e. 

■  *|V2|V2»  (52) 

where  |v2|  is  the  magnitude  of  the  depth  averaged  lower  layer 
current.  A  constant  value  of  2.5  x  10-3  is  assumed  for  the 
coefficient  k.  The  lower  layer  current,  |v2j  is  given  in  terms  of 
the  modal  velocities  by  the  relation 

V2  *  (l/D)Ve  -  (l/H2)Vit  (53) 

where  Vfi  and  are  the  external  and  internal  velocity  vectors, 
respectively. 

The  bottom  stress  for  the  external  mode  computation  at  odd  time 
steps  is  coded  as 

Tb  =  2.5  x  10~3ht { [ ( 1/D ) Ug ( I , J )  -  (1/H2)u”(I,J)]2 

♦  [(l/D)Vg(I, J)  -  (l/H2)v"(i,J)]2}^{-(l/H2)Ui(i,j)}.  (54) 


Upon  substituting  (53)  in  (52)  both  Ue  and  U|  from  the  previous  time 


step  are  used  in  |v2|  while  only  the  previous  value  of  is  used  in 
V2.  This  implies  that  the  coefficient  of  Ue  on  the  left  hand  side  of 
(42)  has  to  be  changed  from  a  constant  value  of  1  to  1  -  2.5  x 
10“3(At/D).  However,  the  tridiagonal  form  of  the  coefficient  matrix 
is  retained. 

At  even  time  steps,  coding  for  the  bottom  stress  is  obtained  by 
directly  interchanging  U  and  V  to  obtain 

Tb  =  2.5  x  10“3At{[(l/D)v"(I,J)  -  (1/H2)V^(I, J)]2 

+  l(l/D)Ug(i,J)  -  (l/H2)u"(i,J) ]2}^{-(1/H2)v"(I, J)} .  (55) 

The  bottom  stress  for  the  internal  mode  computation  is  depicted  as 

Tfa  =  2.5  x  10-3At{[(l/D)Ug(I,J)  -  (1/H2)U^(I, J) ]2 

♦  [(l/D)Vg(I,J)  -  (l/H2)vJ(i,J)]2}^{(l/D)Ug(i,J)}  (56) 

Tb  =  2.5  x  10-3dt{[(l/D)Vg(i,J)  -  (1/H2)v"(I,J)]2 

♦  [(l/D)Ug(i,J)  -  (l/H2)Ui(I,J)]2}^{(l/D)Vg(I,j)}  (57) 

for  odd  and  even  time  steps,  respectively. 

■4 

The  velocity  vector  W  in  the  interface  stress  is  the  internal 
mode  volume  transport  per  unit  width  multiplied  by  the  appropriate 
depth 

W  =  (D/H^H2)Qj,  (58) 

A  constant  value  of  2.0  x  10”3  is  assumed  for  the  friction 
coefficient,  k.  The  finite  difference  form  of  this  stress  is,  for 
odd  time  steps, 
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and 


Tj.  *  2.0  x  10"5{[Ui(I,J)]2+[Vi(l,j)32}*{uJ(I,J)+v"(I,j)}, 

(59 

=  2.0  x  10“5{[vJ(I,J)]2+Cu5(I,J)]2}^{V^a,J)+U^(I,J)}, 
for  even  time  steps, 
d)  The  coupling  terms 

The  coupling  terms  are  coded  based  on  the  premise  that  the 
energy  transferred  between  modes  due  to  coupling  must  be  balanced. 

The  energy  equations  are  formed  by  multiplying  (37a, b)  and  (38a, b)  by 
Qe/D,  g 4/e,  Qj/r^/  and  eg^,  respectively,  and  thus  the  coupling  terms 
appear  as 

-eg^iQev(Hi/D)  (60a; 

-«g^eQiv(H1/D)  (60b: 

*^iQiV(Hl/D)  (60c; 

cg^eQev(H1/D)  (60d: 

For  simplicity,  consider  a  one  dimensional  channel  in  cartesian 
coordinates  as  sketched  in  figure  4.  The  solid  boundary  is  at  point 
1  while  there  are  three  possible  conditions  at  point  7;  i)  a  solid 
boundary,  ii)  an  open  port,  and  iii)  an  interior  point  where  the 
interface  intersects  the  bottom.  Since  the  depth,  and  hence  the 
ratio  Hj/D,  are  defined  at  every  field  variable  location  the 
gradients  [|;^(H1/D)  ]  are  defined  at  midpoints  between  U  and  it 
locations  (refer  to  points  A,  B,  C,  ...  in  Fig.  4).  As  an  example, 
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[|3t<Hi/D)3A  =  (H1/D)?  -  (H1/D)gl 

where  numeral  s'-bscripts  indicate  the  points  where  the  ratios  are 
computed . 

Por  the  external  mode  computations,  the  continuity  and  momentum 
equations  are  alternately  applied  at  points  2  through  6,  while  the 
internal  mode  computations  start  at  point  4  in  the  same  sequence. 

The  summation  of  the  coupling  terms  for  one  complete  operation,  after 
eliminating  like  terms,  is 

-«g[ue]3  JUfiljlfev d>]e} 

♦«9t*i  16  ^nUe]7[|5(H1/D)]A} 

(62) 

-«g[^e^2  +  [Ui]i[|i(H1/D)]F} 

~«g[^e]6  i*[uiM$x(Hi/D)  W  • 

The  internal  mode  computations  start  at  point  4  as  discussed  earlier. 
Hence  the  first  and  third  terms  in  (62)  vanish.  The  residual  terms 
depend  upon  the  condition  at  point  7.  if  point  7  is  a  solid 
boundary,  then  [Ue]7  =  [lh]7  =  0  and  the  net  energy  transferred 
between  modes  is  rero.  The  second  possibility  that  point  7  is  an 
open  boundary  requires  that  either  the  depths  at  point  6  and  7  are 
set  equal  or  the  gradient,  [|x(Hl/D)^A  is  set  t0  zer0*  In  £act  1116 
two  choices  imply  one  another.  It  is  more  rational  to  set  the 
gradient  Hj/D  to  aero  since  all  these  gradients  are  computed  only  for 
the  purpose  of  evaluating  the  coupling  terms.  If  the  interface 
intersects  the  bottom  at  point  7,  then  the  external  mode  computations 


are  continued  until  a  solid  boundary  or  an  open  port  is  encountered, 
while  the  internal  mode  computations  stop  at  point  6.  It  immediately 
follows  that  the  last  term  in  (62)  vanishes.  The  coupling  term 
obtained  by  applying  the  external  momentum  equation  at  point  7  is 

-*g[Ue]7  iC[<ri36C|5<H1/D)lA},  (63) 

which  exactly  balances  the  only  remaining  term  in  (62).  Note  that 
the  second  part  of  the  coupling  term  in  (63)  is  omitted  since  it 
involves  which  is  zero.  It  can  be  shown  from  (59a,b)  that  all 
the  coupling  terms  are  zero  beyond  point  7.  As  a  result,  the  net 
transfer  of  energy  between  modes  is  zero. 

In  summary,  with  this  form  of  coding,  employing  the  average  of 
the  products  U[-|j(H1/D)]  or  ^[^(Hj/D)],  the  only  constraint  needed 
to  fulfill  the  energy  requirement  is  that  the  gradient  of  (H^/D) 
along  the  points  just  inside  the  open  boundary  be  zero. 

Therefore,  the  coupling  terms,  x^»  and  in  the  implicit 
computations  of  the  external  mode  momentum  and  mass  conservation 
equations,  (41)  and  (43),  are  coded  as  follows: 

xi  =  f7x  ,J)[(H1/D(I,J))-(H1/D(I-JfJ))] 

'*-^(I'J4)^Hl/D(I4'J))"(Hl/D(I»J)>^'  (64a) 

*i  =  J«{ux{uJ(I-i.J)[<H1/D(I#J)>-(H1/D(I-J#J))] 

J)  [(Hi/DCI^,  J>  >-<H1/D(I,J)  >  3) 
♦uy{V^(I,j4)«CJ^)[(H1/D(I,j4))-(H1/D(I,J))] 
♦Vj(I.J-i)«(J-J)[{Hl/D(I,J))-(Bl/D(I,J-i))])} .  (64b) 

For  the  explicit  calculations  at  odd  time  steps,  x^,  in  (46)  becomes 
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+^(I-J,j)5(j-J)[(H1/D(l,j))-(H1/D(lrj-J))]}.  (65] 

Note  that  the  values  of  all  the  field  variables  at  previous  time 
steps  are  used  in  the  expressions  of  the  coupling  terms.  This 
implies  that  the  codes  for  these  terms  at  even  time  steps  are  the 
same  as  those  at  odd  time  steps  except  for  the  sequence  in  which  they 
are  applied.  Eqs.  (65),  (64b),  and  (64a)  are  the  codes  employed  for 
implicit  computations  of  V  and  &  and  for  explicit  computations  of  U, 
respectively. 

Sequential  coding  of  coupling  terms  for  the  implicit  and 
explicit  internal  mode  computations  at  odd  time  steps,  are 

xe  =  -g[«h1H2/D(l,j)]ukJ{*J(i-J,j)[(H1/D(I,j))-(H1/D(l-JtJ))] 

+^(I,J-4)[(H1/D(l4»J))-(H1/D(l,J))]},  (66a) 

le  =  4^x^ue(I_i'J>^Hl/D(I»J)>"(Hl/D<I_i»J))^ 

■*-U^(l4#J)[(H1/D(I^,J))-(H1/D(I,J))]} 

♦uy{V^(I,J-^)d(J^)[(H1/D(I,j4))-(H1/D(I,J))] 

♦Vgd,  J-^)0(J-^)  [(H^/D(I, J))-(Hj/D(I,J-^)) ]}} ,  and  (66b) 

xe  *  -g[«H1H2/D(i,j)3uyJ{^(i,j4)«(j4)f(Hi/D(I'J4))) 

-(H1/D(I,J))]>^(I-^,J)$(J-^)  [(Ni/Dd.J)  )-(H1/D(I,J-^)  )]}  .  (67) 

The  same  considerations  for  the  coding  at  even  time  steps  as 
discussed  above  also  apply. 

It  is  noteworthy  that  the  numerical  coding  of  the  coupling  terms 
employed  is  the  only  possible  form,  on  the  basis  of  energy  transfer, 


for  the  spatial-stagger ed  grid  system  used  in  this  study. 

Furthermore,  the  energy  transfer  is  balanced  globally,  but  not 
locally.  In  other  words,  a  balance  is  obtained  when  considering  the 
entire  domain  but  not  at  any  individual  grid  block. 

e)  Initial  and  boundary  conditions 

The  model  was  taken  initially  at  rest.  The  initial  positions  of 
the  storms  were  in  the  Cayman  Sea  for  the  hurricanes  of  record  and 
for  all  except  one  synthetic  storm  employed  in  the  parametric  study 
which  will  be  discussed  in  Chapter  IV. 

Specified  volume  transport  or  height  anomalies  are  employed  as 
boundary  conditions.  Along  the  solid  coastal  boundaries  the  normal 
component  of  volume  transport  is  specified  as  zero.  At  the  open  grid 
elements  water  levels  are  placed  in  equilibrium  with  the  inverted 
barometric  pressure  for  the  external  mode  and  zero  for  the  internal 
mode. 


3.  Wind  and  pressure  forcing 
a)  Analytical  wind  model 

Fracticaliy,  there  are  two  methods  of  portraying  hurricane  wind 
fields  on  the  computing  grid.  One  method  is  to  digitize  the  surface 
charts  that  are  available.  The  data  are  sampled  at  stipulated  time 
intervals  and  grid  points  and  then  interpolated  in  space  and  time  to 
provide  the  necessary  information  to  the  model.  This  method  requires 
detailed  surface  charts  of  wind  throughout  the  simulation  period.  In 
addition,  this  is  a  laborious  technique  to  apply,  especially  for  this 


study  which  requires  these  data  over  long  time  periods.  The 
alternative,  which  is  used  in  this  study,  is  to  derive  the  required 
forcing  fields  from  a  parametric  analytical  model. 

Schwerdt  et  al.  (1979)  developed  a  model  for  surface  wind  fields 
associated  with  hurricanes  for  the  Gulf  and  East  Coasts  of  the  United 
States  which  is  commonly  known  as  National  Weather  Service  Model, 
NWS-23.  This  model  was  employed  to  reconstruct  radial  wind  profiles 
for  hurricane  Carla.  Comparisons  between  the  model  results  and  the 
observed  profiles  revealed  that  this  model  does  not  satisfactorily 
depict  the  observed  Carla  winds  in  the  far  field.  It  is  remarked 
that  effect  of  the  storm  forward  speed  is  ignored 

In  NWS-23  there  are  two  variables,  the  maximum  wind,  Vm,  and  the 
radius  to  maximum  wind,  R,  that  determine  the  wind  speed  for  a 
stationary  storm.  Therefore,  in  order  to  obtain  good  comparison  of 
the  winds  in  the  far  field,  either  Vffl  or  R,  or  both  have  to  be 
changed  and  this  inevitably  deteriorates  the  winds  near  the  center  of 
the  storm. 


Holland  (1980)  proposed  a  new  model  that  has  two  parameters,  C 
and  k  that  independently  define  the  location  of  the  maximum  winds  and 
the  shape  of  the  wind  profile,  respectively.  The  model  thus  allows 
the  adjustment  of  wind  speed  in  the  far  field  through  the  parameter 
k,  without  changing  R.  The  gradient  wind  profile  is  given  as 


[Ck(Pn  -  Pc)exp(-C/r!t)/parJt  ♦  (£rf)2]-  -  £rf. 


where  wg  is  the  gradient  wind  at  a  distant  r  from  the  center,  f  is 
the  Coriolis  parameter,  pa  is  the  the  air  density  (assumed  constant). 


Pc  is  the  central  pressure,  and  Pn  is  the  ambient  pressure. 


Holland  determined  C  and  k  by  fitting  the  pressure  profile 

P  =  Pc  +  <Pn  -  Pc)exp(-C/rk).  (69 


However,  he  pointed  out  that  this  approach  would  underestimate  the 
peak  winds  and  that  wind  observations,  if  available,  should  be  used 
directly. 

For  hurricane  Carla,  there  were  detailed  surface  charts 
available.  Equation  (68)  was  modified  for  the  direct  approach  as 
follows:  First,  C  was  eliminated  by  considering  that  the 
cyclostrophic  wind 


Wc  =  [Ck(Pn  -  Pc)exp(-C/rk)/park] 


J 


(70 


was  a  more  appropriate  representation  of  the  maximum  wind.  Upon 
taking  dwc/dr  and  setting  it  equal  to  zero  at  r=R,  we  obtain  the 
relation 


R  =  C1//k  or  C  =  Rk 


(71 


Since  detailed  surface  charts  allow  a  good  estimate  o'  the  maxi mum 
wind  and  radius  of  maximum  wind,  the  pressure  drop  (Pn~Pc)  was 
represented  in  terms  of  the  maximum  wind  by  substituting  (71)  in  (70) 
to  obtain 


W 


m 


[(k/pae) (Pn 


(72 


Using  (71)  and  (72),  (68)  reduces  to 

Wg  =  [(R/r)kW^cxp(l-(R/r)K)  ♦(^rf)<i]*:  -  irrf. 


<71 


It  should  be  remarked  that  for  k=l,  this  wind  model  is 


equivalent  to  the  NWS- 23  model.  Another  important  point  to  note  is 
that  the  azimuthal  variation  of  k  automatically  results  in  an 
asymmetry  of  the  wind  field  due  to  the  translation  of  the  storm. 

Prom  the  NWS-23  model,  this  asymmetry  has  to  be  taken  into  account  by 
augmenting  the  maximum  wind  speed  with  the  forward  speed  of  the  storm 
and  the  cosine  of  the  angle  that  depends  upon  the  location  of  the 
maximum  wind  relative  to  the  storm  path. 


b)  Analytical  pressure  model 

The  pressure  profiles  obtained  from  (69)  using  k  obtained  from 
fitting  the  velocity  profiles  were  not  in  good  agreement  with  the 
observed  pressure  profiles.  This  confirms  Holland's  remark  as 
mentioned  earlier.  Therefore,  the  pressure  profile 

p  =  Pc  ♦  (Pn-Pc)e(“R/r),  (74 

which  is  equivalent  to  (69)  with  k=l,  was  employed  in  this  study. 


CHAPTER  III 


SIMULATIONS  OF  HURRICANES  OF  RECORD 

1.  Selected  hurricanes  of  record 

Eligible  hurricanes  of  record  for  the  purposes  of  this  study 
were  considered  as  those  storms  for  which  the  historical  information 
required  by  the  numerical  model  was  available.  This  includes  the 
time  series  of  the  central  pressure,  radius  of  maximum  wind  and 
surface  charts  of  wind  and  pressure.  These  charts  should  cover  the 
entire  Gulf  of  Mexico  and  Cayman  Sea  from  the  time  at  which  the  storm 
center  was  outside  the  Gulf  to  sometime  after  landfall.  Additional 
important  data  required  includes  the  water  level  associated  with 
these  storms  at  stations  around  the  Gulf,  in  both  the  United  States 
and  Mexico. 

Of  the  26  hurricanes  spanning  the  period  1950-1980  which  were 
examined  as  potential  hurricanes  of  record,  only  hurricanes  Carla  in 
1961  and  Allen  in  1980  had  sufficient  observations  of  the  forcing 
fields  and  the  response  histories.  These  are  the  hurricanes  of 
record  used  for  verification  purposes  in  this  study. 

Hurricane  Carla  was  an  exceptionally  slow  moving  storm  with  an 
average  forward  speed  of  13  km/h.  Figure  5  shows  the  path  indicated 
by  serial  positions  of  the  storm's  center  at  six  hour  intervals. 

Carla  reached  hurricane  stage  at  1200  GMT  6  September  1961  and  th© 
center  entered  the  Gulf  through  Yucatan  Strait  at  approximately  1500 
GMT  7  September  with  a  central  pressure  of  970  mb.  As  it  moved  to 
the  northwest,  it  continously  deepened  and  reached  a  minimum  central 


Track  of  hurricane  Carla  from  1200  GMT  6  September  to  1800  GMT  11  September  1961 


pressure  of  935  mb  at  1200  GMT  11  September.  The  time  sequence  of 
the  central  pressure  from  September  4  to  September  13,  Fig.  6,  was 
presented  by  Dunn  et  al.  (1962).  During  the  period  1200  GMT  9 
September  through  1800  GMT  12  September,  the  average  radius  of 
maximum  wind  was  approximately  40  km  with  a  slight  increase  observed. 
Upon  entering  the  Gulf,  Carla  had  winds  of  60  km/h.  As  it  drifted 
northwest,  then  west  northwest  it  increased  in  size  with  cyclonic 
winds  observed  over  the  entire  Gulf.  Maximum  wind  speeds  of  82  km/h 
were  observed  inland  as  Carla  approached  the  coast.  Carla  made 
landfall  near  Pass  Cavallo  at  2100  GMT  11  September  and  by  12 
September  it  was  positioned  north  of  Waco.  Hurricane  Carla  has  the 
distinction  of  being  the  best  documented  storm  in  history. 

Hurricane  Allen  was  a  fast  moving  storm  compared  to  Carla.  Its 
average  forward  speed  in  the  Cayman  Sea  was  35  km/h.  Its  forward 
speed  decelerated  to  about  30  km/h  as  it  moved  west-northwest  across 
the  Gulf.  The  center  of  the  storm  entered  Yucatan  Strait  at 
approximately  1800  GMT  7  August,  1980.  The  hurricane  center  crossed 
the  coastline  near  Brownsville  at  0700  GMT  10  August.  Figure  7  shows 
the  path  of  Allen  from  the  Cayman  Sea  until  landfall. 

The  time  sequence  of  Allen's  central  pressure,  as  shown  in  Fig. 
8,  was  obtained  from  Lawrence  and  Pelissier  (1981).  There  were  three 
cycles  of  a  50  mb  fluctuation  in  the  six-day  period  from  4-10  August. 
The  first  minimum  central  pressure  of  911  mb  occurred  early  on  5 
August  when  the  storm  was  approximately  370  km  south  of  Puerto  Rico. 
The  next  minimum  of  899  mb  was  measured  at  1742  GMT  7  August  when  the 
eye  just  passed  through  the  strait.  Allen  deepened  for  the  third 
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ig.  6.  Time  sequence  of  Carla  central  pressure  4 
1961  (after  Dunn  et  al.,  1962). 


-  13  September 


time  when  it  moved  toward  Texas  coast  and  reached  a  909  mb  low  early 
on  9  August  local  time. 

2.  Meteorological  data 

The  Hydrometeorological  Section  of  the  U.S.  Weather  Bureau 
provided  surface  charts  of  hurricane  Carla  winds  and  pressure  for  the 
entire  Gulf  of  Mexico.  These  charts  were  available  at  6  h  intervals 
from  1200  GMT  9  September  to  1200  GMT  10  September  and  at  3  h 
intervals  thereafter .  Prior  to  this  period,  the  surface  pressure 
charts  that  covered  North  America,  obtained  from  the  National 
Climatic  Center,  provided  observed  surface  winds  from  land  stations, 
a  bouy  and  ships  of  opportunity. 

Ho  and  Miller  (1980)  presented  several  surface  wind  charts  for 
the  period  when  Allen  was  in  the  western  Gulf.  The  coverage  of  these 
charts  is  limited  to  the  western  Gulf  only.  Despite  the  lack  of 
surface  wind  charts  however,  the  time  sequence  of  the  central 
pressure  and  the  available  surface  pressure  charts  provided  enough 
information  to  construct  surface  wind  fields.  However,  this 
information  was  not  sufficient  to  provide  the  same  degree  of  detail 
as  obtained  with  hurricane  Carla. 

a)  Surface  wind  for  hurricane  Carla 

The  hurricane  Carla  surface  wind  charts  were  analyzed  to 
determine  the  location  of  the  eye  and  R,  wffi  and  k  which  were  required 
for  the  computation  of  wind  fields. 
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The  position  of  the  eye  and  the  radius  of  maximum  wind  were 
first  determined.  To  account  for  the  asymmetry  of  the  wind  fields,  k 
was  determined  by  fitting  the  observed  radial  wind  profiles  along 
different  sectants  around  the  eye.  A  total  of  8  profiles  was  fitted 
for  each  chart.  Figure  9  shows  sections  where  these  profiles  were 
fitted.  In  each  case,  the  profile  was  plotted  by  digitizing  the 
radial  distance  from  the  center  to  the  isovels  on  the  observed  chart. 
The  maximum  wind  speed,  Wm,  was  determined  from  the  plot  at  the 
distance  R  from  the  center.  The  radial  wind  profile  was  computed 
from  (73)  by  substituting  R,  Wm  and  assuming  k=l.  The  resulting 
profile  was  compared  to  the  observed  and  k  was  adjusted  to  obtain  the 
best  possible  agreement. 

Prior  to  1200  GMT  9  September,  during  which  there  were  no 
regional  Carla  surface  wind  charts,  the  analysis  to  obtain  R,  Wm  and 
k  depended  upon  the  surface  pressure  charts  from  the  National 
Climatic  Center.  Eq.  68  was  rewritten,  by  using  (71),  as 

W  =  [(R/r)k(k/pa)(Pn-Pc)e"(R/r)k  +  (£rf)2]^  -  Jrf.  (75 

The  time  sequence  of  the  radius  of  maximum  wind  was  extrapolated  back 
in  time  using  its  relation  to  the  central  pressure  which  is 
summarized  in  NWS-23.  The  radial  distances  to  those  points  where 
observed  winds  were  reported,  were  digitized  from  the  surface  chart. 
The  far  field  pressure  was  determined  from  the  first  cyclonically 
curved  isobar.  Assuming  pa  =  1.15  Kg/m2  and  k»l,  the  wind  speed 
computed  from  (75)  was  compared  to  the  observed  value  and  again  k  was 
adjusted  to  obtain  minimum  error.  Once  k  was  determined,  was 


computed  from  (72). 

Attempts  were  made  to  obtain  a  simple  expression  for  both  k  and 
Wm  as  a  function  of  the  azimuthal  angle  (measured  from  the  north) 
such  that  these  two  parameters  could  be  internally  computed. 

However,  no  simple  relationship  could  be  established  due  to  the 
complexities  of  the  patterns  of  azimuthal  variation  of  k  and  Wm 
which,  in  addition,  varied  irregularly  for  each  map  time.  Therefore, 
these  parameters  have  to  be  specified  serially.  Furthermore,  linear 
interpolation  in  time  to  obtain  these  parameters  for  intermediate 
time  step  computations  constrained  the  values  to  the  same  set  of 
azimuths  (relative  to  heading)  throughout  the  simulation  period.  As 
a  result,  an  a  priori  linear  interpolation  in  space  was  employed  to 
get  k  and  Wm  as  continuous  functions  of  azimuth.  A  standard  set  of 
azimuths  was  then  selected. 

Finally,  the  surface  wind  fields  for  each  map  time  were 
constructed  using  the  interpolational  routine  employed  in  the  model. 
Comparison  between  the  computed  and  observed  wind  fields  was  made  and 
final  adjustments  of  k,  if  necessary,  were  decided.  Figures  10  and 
11  are  examples  of  the  computed  wind  fields  at  0000  GMT  7  September 
and  1200  GMT  9  September,  respectively.  The  solid  circles  indicate 
the  observed  wind  speed  from  the  surface  charts. 

b)  Surface  wind  for  hurricane  Allen 

The  inadequate  coverage  in  both  space  and  time  of  surface  wind 
charts  for  hurricane  Allen  aborted  attempts  to  fully  analyze  these 
charts  as  practiced  in  hurricane  Carla.  Therefore,  NWS-23  was  used 
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Hurricane  Carla  surface  wind  (in/s)  at  0000  GMT  7  September  1961,  as  obtained  from  Holland’ 
tie  solid  circles  indicate  the  observed  wind  speed  from  the  surface  chart. 


mmMmmMmummmm 


Hurricane  <larla  surface  wind  (m/s)  at  1200  GMT  9  September  1961,  as  obtained  from  Holland* 
he  solid  circles  Indicate  the  observed  wind  speed  from  the  the  surface  chart. 


to  determine  R  and  Wm  from  Pc  which  were  required  to  construct  the 
wind  fields.  The  graphical  relationship  between  the  radius  of 
maximum  wind  and  the  central  pressure  was  analyzed  using  simple 
linear  regression  to  establish  a  functional  relationship.  The  radius 
of  maximum  wind  was  then  determined.  The  maximum  wind  speed  for  a 
stationary  storm  Wms,  is  given  as 

Wms  =  0.9(l/pe)z(Pn-Pc)2  -  (Rf/2)  (76) 

Since  there  were  no  surface  charts  to  analyze  for  the  parameter  k,  a 
constant  value  of  1  was  assumed.  As  discussed  earlier,  the  asymmetry 
of  the  wind  field  was  achieved  by  augmenting  the  maximum  wind  speed 
by 

Wm  =  Wms  +  1.5(Vf°-63)(T§-37)coso,  (77) 

where  o  was  the  angle  between  track  direction  and  the  surface  wind 
vector  and  T0  =  0.514791  for  wind  speed  in  m/s.  The  track  of 
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hurricane  Allen  was  approximately  285  (cf.  Fig.  7,  p.  41)  relative 
to  true  north.  Assuming  that  the  maximum  wind  occurred  at  right 
angles  to  the  right  of  the  track,  then  o  was  zero  at  25  relative  to 
true  north.  Eq.  (78)  was  used  to  compute  the  maximum  wind  speed  for 
each  grid  point  (I,J)  and  thus  the  surface  wind  speed  at  all 
computational  points  can  then  be  evaluated. 


3.  Tide  gauge  data 


A  total  of  9  and  13  tide  gauge  stations  for  hurricane  Carla  and 
Allen,  respectively,  were  chosen  to  provide  the  observed  response  in 
the  Gulf  of  Mexico.  Figure  12  shows  locations  of  these  tide  gauge 
stations.  The  Tides  Prediction  Branch  of  the  National  Oceanic  and 
Atmospheric  Administration  (NOAA)  and  the  National  Autonomous 
University  of  Mexico  (UNAM)  provided  hourly  water  level  for  a  period 
of  several  weeks  before  and  after  the  hurricanes.  Filtered  versions 
(using  a  40  h  lowpass  filter)  ar<e»  shown  in  Figs.  13-15  for  some 
selected  series  obtained  during  hurricane  Carla.  These  filtered 
plots  essentially  remove  the  tides.  The  arrows  in  the  figures 
indicate  the  time  at  which  the  center  of  Carla  entered  the  Gulf 
through  Yucatan  Strait.  There  were  indications  of  a  gradual  rise  of 
water  level  well  before  the  peak  surge,  a  possibility  of  forerunner 
surges  at  all  9  stations.  Note  also  the  degree  of  background 
variation  unrelated  to  the  hurricane  several  weeks  prior  to  the 
hurricane. 

The  same  analysis  as  applied  to  recorded  water  level  for 
hurricane  Carla  were  used  for  the  records  obtained  obtained  during 
hurricane  Allen.  Figures  16-18  show  the  filtered  data  at  some 
selected  stations.  The  presence  of  an  initial  rise  in  water  level  is 
again  observed  at  all  stations. 

The  lowpass  filtering  illustrated  above  is  known  to  smooth  out 
and  broaden  the  peak  surge.  To  properly  r amove  tidal  signals  from 
recorded  data  by  the  harmonic  method  requires  a  suitable  number  of 
constituents  to  insure  proper  phasing.  Unfortunately,  only  the  Key 


ig.  13.  Observed  water  levels  at  Galveston  during  hurricane  Carla,  1961.  The  arrow  indicates  the 
time  at  which  Carla  entered  the  Gulf  through  Yucatan  Strait.  The  datum  is  gauge  mean  low  water. 


Observed  water  levels  at  Port  Isabel  during  hurricane  Carla,  1961.  The  arrow  indicates  the 
which  Carla  entered  the  Gulf  through  Yucatan  Strait.  The  datum  is  gauge  mean  low  water. 


Fig.  15.  Observed  water  levels  et  Key  West  during  hurricane  Carla,  1961.  The  arrow  indicates  the 
time  at  which  Cat  la  entered  the  Gulf  through  Yucatan  Strait.  The  datum  is  gauge  mean  low  water. 


6.  Observed  water  levels  at  Galveston  during  hurricane  Allen,  1580.  The  arrow  indicates  the 
at  which  Allen  entered  the  Gulf  through  Yucatan  Strait.  The  datum  is  gauge  mean  low  water. 


Fig.  17.  Observed  water  levels  at  Port  Isabel  during  hurricane  Allen,  1980.  The  arrow  indicates  the 
time  at  which  Allen  entered  the  Gulf  through  Yucatan  Strait.  The  datum  is  gauge  mean  low  water. 


ed  water  levels  at  Key  West  during  hurricane  Allen,  1980.  The  arrow  Indicates  the 
lien  entered  the  Gulf  through  Yucatan  Strait.  The  datum  is  gauge  mean  low  water. 


West  and  Galveston  stations  have  enough  information  for  the  actual 
computation  of  the  tide.  Therefore,  at  each  station,  only  the  mean 
water  level  during  the  period  considered,  excluding  the  peak  surge, 
was  removed  from  the  raw  data.  Thus  the  recorded  hydrographs  which 
are  compared  to  those  computed  contain  the  tidal  signal,  but  its 
amplitude  is  small  and  does  not  unduly  complicate  the  comparison. 

4.  Simulation  procedure 

The  initial  model  wind  and  pressure  fields  were  increased  in 
amplitude  linearly  from  zero  to  their  actual  initial  values  over  a  48 
hour  prototype  time  period.  During  this  period  Carla  was  kept 
stationary  at  19.04°N  and  85.15°W  which  was  its  position  at  1200  GMT 
6  September.  The  simulation  was  carried  out  to  0000  GMT  13 
September,  approximately  28  h  after  landfall. 

The  simulation  for  hurricane  Allen  starts  at  0000  GMT  7  August 
and  ends  at  0000  GMT  12  August.  The  initial  position  of  the  eye  of 
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hurricane  Allen  during  start  up  was  at  20.10  N  and  81.90  W. 

At  each  time  step  the  east  and  north  component  of  wind  stress 
and  the  atmospheric  pressure  (the  inverted  barometric  height)  were 
computed  at  each  grid  point.  The  sequence  of  wind  and  pressure 
computations  consisted  of  two  linear  interpolations.  First,  the 
position  of  the  eye  and  the  other  parameters  at  each  time  step  were 
linearly  interpolated  from  two  appropriate  sets  of  values  (6  h 
apart).  The  radial  distance  from  the  eye  to  each  individual  grid 
point  (I,J)  was  computed.  The  inverted  barometric  height  was  then 
computed  from  (74).  The  angle  between  the  line  joining  the  eye  to 
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the  grid  point  (I,J)  and  true  north,  8,  was  determined.  The  pair  of 
azimuthal  angles  that  embraced  8  was  sorted.  The  values  of  k  and  Wm 
associated  with  the  two  angles  so  determined  were  used  in  the  linear 

interpolation  to  obtain  their  values  at  grid  point  (I,J). 

0 

A  constant  inflow  angle  of  20  was  assumed  in  all  simulations 
except  Carla  for  the  decomposition  of  the  wind  speed  before  computing 
the  wind  stress  components.  The  computed  water  level  at  Galveston 
during  the  early  stage  of  Carla  simulation  using  constant  inflow 
angle  was  lower  than  the  observed.  Careful  examination  of  the 
surface  wind  charts  at  1200  GMT  9  September  revealed  that  there  was  a 
region  along  the  Texas  shelf  where  the  cross  shelf  wind  reversed 
direction  indicating  a  negative  inflow  angle.  An  example  of  an 
analyzed  wind  map  given  by  Miyasaki  (1963),  also  shows  negative 
inflow  angles  along  the  Texas  shelf.  Therefore,  the  inflow  angle  was 
allowed  to  vary  as  a  function  of  the  radial  distance  using  the 
empirically  determined  formula 
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for  r  <  R 


1' 


(78) 


02-(r-R1)*O.5*A/r  for  r  >  R^, 


where  A  is  a  maximum  inflow  angle,  20*,  R^  is  the  radial  distance  at 
which  the  inflow  angle  was  zero,  and 


0:  *  [(R1/r)A2exp(l-(R1/r))*(Jrf )2J  -  jrf, 
ej  •  [(Rj/r)*  A2expU-(R1/r)K)*($r£)2)  -  Jrf, 
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K  =  r/50  . 


Figure  19  shows  the  adopted  inflow  angle  profile.  This  radially 
varying  inflow  angle  was  applied  only  from  the  begining  of  the 
simulation  to  0000  GMT  10  September.  A  constant  inflow  angle  was 
resumed  after  this  period  since  there  was  no  other  evidence  of 
negative  inflow  angles  on  the  remaining  available  surface  charts. 

Results  of  the  simulation  were  sampled  every  24  hour  prototype 
time.  These  included  digital  fields  and  map  plots  of  the  height 
anomalies  and  currents  for  both  modes,  and  surface  currents  whose 
values  were  retrieved  from  the  modal  currents.  At  the  end  of  the 
simulation,  the  computed  and  observed  hydrographs  at  the  tide 
stations  were  plotted.  Contours  of  the  peak  surge  on  the  continental 
shelf  for  the  entire  Gulf  of  Mexico  and  for  the  Texas-Louisiana  shelf 
were  also  plotted. 

5.  Results  of  Carla  simulation 

Figures  20  through  28  show  the  computed  (solid)  and  recorded 
(dashed)  hydrographs  at  the  stations  used  in  this  study.  The  overall 
comparisons  are  fair,  especially  during  the  first  half  of  the 
simulation  period.  There  is  not  ouch  activity  in  the  southwestern 
portion  of  the  Gulf  (Campeche  Say)  as  revealed  from  the  hydrographs 
from  the  three  Mexican  stations  (Figs.  26  through  28).  The  simulated 
hydrograph  at  Key  West  has  a  slow  oscillation  with  a  period  of 
approximately  7  days.  The  mean  of  this  long  period  oscillation  lies 
close  to  the  mean  of  the  recorded  water  ievel.  This  signax  is  also 
observable  in  the  St.  Petersburg  computed  hydrograph.  However,  the 
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Fig*  20.  Computed  (solid)  and  observed  (dashed)  hydrographs  at  Key  West  during  hurricane  Carla 
6-13  September  1961.  The  datum  is  mean  sea  level. 


22.  Computed  (solid)  and  observed  (dashed)  hydrographs  at  Pensacola  during  hurricane  Carla 
13  September  1961.  The  datum  is  mean  sea  level. 
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Fig.  23.  Computed  (solid)  and  observed  (dashed)  hydrographs  at  Grand  Isle  during  hurricane  Carla 
6-13  September  1961.  The  datum  is  mean  sea  level. 
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24.  Computed  (solid)  and  observed  (dashed)  hydrographs  at  Galveston  during  hurricane  Car 
13  September  1961.  The  datum  is  mean  sea  level. 
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13  September  1961.  The  datum  is  mean  sea  level. 
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Computed  (solid)  and  observed  (dashed)  hydrographs  at  Veracruz  during  hurricane  Carla 
September  1961.  The  datum  is  mean  sea  level. 
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Computed  (solid)  and  observed  (dashed)  hydrographs  at  Carmen  during  hurricane  Carla 
;ptember  1961.  The  datum  is  mean  sea  level. 
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Computed  (solid)  and  observed  (dashed)  hydrographs  at  Progreso  during  hurricane  Carla 
September  1961.  The  datum  is  mean  sea  level. 
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computed  water  level  at  St.  Petersburg  is  higher  than  the  recorded 
throughout  the  period.  This  m:.ght  be  due  to  the  fact  that  this  tide 
gauge  station  is  located  inland  and  not  on  the  open  coast.  The 
departure  of  the  two  hydrographs  at  Grand  Isle  is  probably  the  result 
of  resuming  the  20°  inflow  angle  which  shifted  the  wind  direction 
from  thu  alongshore  or  slightly  onshore  to  slightly  offshore.  The 
comparison  for  Galveston  is  better  than  that  obtained  by  Miyasaki 
(1963),  presumably  because  of  the  special  care  to  match  the  observed 
wind  fields  for  the  northwest  Gulf  Coast.  At  Port  Isabel,  the 
discrepancy,  which  resembles  that  at  Grand  Isle  but  with  a  larger 
deviation,  is  possibly  due  to  the  difference  in  locations  of  the  tide 
gauge  station  and  the  point  where  the  computed  water  level  was 
sampled.  The  actual  site  of  the  tide  gauge  station  is  inside  a  semi- 
enclosed  embayment,  as  shown  in  Fig.  29.  As  a  result,  the  effect  of 
the  offshore  component  of  the  wind,  which  prevailed  at  this  station 
from  approximately  0000  GMT  10  September  onward,  is  limited  while  the 
southerly  component  produces  a  set  up  within  the  constricted  lagoon 
at  the  gauge  site.  On  the  contrary,  the  computed  hydrograpn,  sampled 
at  half  the  grid  size  away  from  the  digitized  coastline,  is  subjected 
fully  to  the  wind  draw-down.  Miyasaki  (1963)  obtained  a  similar 
result  at  this  station. 

Despite  these  discrepancies,  the  primary  interest  is  the 
comparison  at  Galveston  where  the  largest  surge,  among  all  sampled 
stations,  occurs.  It  should  be  pointed  out  that  the  maximum  computed 
water  level  is  not  at  Galveston  but  at  the  grid  point  (8,42),  located 
one  grid  block  northeast  of  the  path  of  Carla  at  the  coast  close  to 
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Fig.  29.  Location  of  Port  Isabel  tide  gauge  station  which  is 
sheltered  by  south  Padre  Island. 


the  location  of  Matagorda  Bay  where  the  highest  high  water  marks  were 
observed  for  this  hurricane.  The  agreement  is  good  over  the  entire 
simulation  period  at  Galveston.  The  effect  of  the  reversing  inflow 
angle  in  the  simulation  is  clearly  pronounced  as  the  computed  water 
level  suddenly  drops  at  approimately  2200  GMT  9  September.  This 
result  is  somewhat  out  of  phase  with  that  indicated  by  the 
observations . 

Pigures  30-32  show  instantaneous  fields  of  the  barotropi.c  and 
baroclinic  height  anomalies  and  surface  currents  at  1200  GWI  10 
September.  Contours  of  the  barotropic  height  (Fig.  30)  clearly 
demonstrate  the  inverted  barometric  effect  around  the  storm  center. 
Upwelling  along  the  hurricane  path,  a  feature  first  investigated  by 
Leipper  (1967)  in  the  wake  of  hurricane  Hilda  in  the  Gulf  of  Mexico 
in  1964,  is  noticeable  in  the  baroclinic  height  field  shown  in  Fig. 
31.  It  should  be  remembered  that  the  negative  contours  of  the 
baroclinic  height  anomalies  correspond  to  upwelling  of  the  interface. 
Geisler  (1970),  in  his  linear  analytic  model,  found  that  the 
baroclinic  response  to  a  moving  hurricane  consisted  of  both  upwelling 
and  inertio-gravity  waves  in  the  lee  of  the  storm.  Figures  31,32 
strikingly  portray  this  wake  oscillation.  Chang  and  Anthes  (1978) 
carried  out  numerical  experiments  to  investigate  the  character  of 
this  wake  produced  by  various  types  of  atmospheric  forcing.  They 
found  that  the  wavelength  is  longer  for  a  faster  moving  storm  and, 
for  the  same  forward  speed  of  the  storm,  is  shorter  at  higher 
lattitude.  In  addition,  within  the  limit  of  SO  h  and  1200  km  time 
and  space  scales  employed,  the  0-effect  (variation  of  f)  did  not 
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Fig.  30.  Computed  barotropic  height  anomaly  field  (meters)  for  hurricane  Carla  at  1200 
September  1961. 


Computed  baroclinic  height  anomaly  field  (meters)  for  hurricane  Carla  at  1200  GMT  10 


hurricane  Carla  at  1200  GMT  10  September  1961 


alter  the  results.  For  an  asymmetric  storm  moving  at  18  km/h,  which 
is  comparable  to  the  translation  speed  of  Carla,  in  an  idealized  f- 
plane  basin,  the  wake  oscillation  has  a  wave  length  of  approximately 
420  km.  The  average  wave  length  obtained  from  this  Carla  simulation 
is  approimately  360  km. 

6.  Results  of  Allen  simulation 

The  lesults  of  the  Allen  simulation  show  a  synchronous 
oscillatory  signal  with  a  period  of  approximately  28  h  in  all  of  the 
simulated  hydrographs.  The  amplitude  of  this  oscillation  is 
approximately  20  cm.  The  existence  of  this  in  phase  signal  in  all 
stations  around  the  Gulf  is  indicative  of  a  Gulf-wide,  Helmholtz  mode 
superimposed  on  other  forced  modes  -  Consequently,  an  ad-hoc  code  to 
compute  the  average  water  level  in  the  Gulf  proper  (designated  as  t?g) 
at  each  time  step  was  added.  Figure  33  shows  the  time  series  of  t)q. 
The  solid  line  represents  tjg  computed  by  averaging  water  levels  at 
each  grid  point  in  the  Gulf,  and  the  dashed  line  is  r?G  computed  from 
the  continuity  equation  using  the  difference  of  volume  transports  at 
the  Florida  and  Yucatan  ports.  The  two  curves  are  practically 
coincident  and  very  well  matched  with  the  signal  present  in  the 
individually  computed  hydrographs  mentioned  above.  As  it  can  be  seen 
in  Fig.  33,  this  mode  persisted  throughout  the  entire  simulation 
period. 

The  agreement  between  the  simulated  and  observed  hydrographs 
(Figs.  34-46)  is  comparable  to  that  obtained  in  the  Carla  simulation. 
This  is  surprising  since  less  information  regarding  the  atmospheric 
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Computed  (solid)  and  observed  (dashed)  hydrographs  at  Key  West  during  hurricane  Allen 
August  1980.  The  datum  is  mean  sea  level. 
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Computed  (solid)  and  observed  (dashed)  hydrographs  at  Cedar  Key  during  harricam 
iqust  1980.  The  datum  is  mean  sea  level. 
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Fig.  39.  Computed  (solid)  and  observed  (dashed)  hydrographs  at  Pensacola  during  hurricane  Allen 
7-12  August  1980.  The  datum  is  mean  sea  level. 
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pig.  41.  Computed  (solid)  artd  observed  (dashed)  hydrographs  at 
7-12  August  1980.  The  datura  Is  mean  sea  level. 
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2  August  I960.  The  datura  is  mean  sea  level 


s.o 


Computed  (solid)  and  observed  (dashed)  hydrographs  at  Veracruz  during  hurricane  Allen 
August  I960.  The  datum  is  mean  sea  level. 


Computed  (solid)  and  observed  (dashed)  hydrographs  at  Carmen  during  hurricane 
August  1980.  The  datum  is  mean  sea  level. 


12  August  1980.  The  datum  is  mean  sea  level. 


forcing  is  available  for  hurricane  Allen  than  for  hurricane  Carla. 

Hydrographs  from  the  Mexican  stations  (Figs.  44-46,  p.  89  -  91) 
reveal  almost  no  response  except  for  tjq  for  the  entire  period 
considered.  The  computed  water  level  at  Progreso  shows  a  response  to 
the  barometric  pressure  drop  during  the  passage  of  the  center  of 
Allen  between  1800  GMT  7  August  to  0006  GMT  8  August.  However,  the 
observed  hydrographs  have  no  indication  of  this  perturbation.  This 
circumstance  also  occurs  in  the  simulation  of  hurricane  Carla  (cf. 
Fig.  28,  p.  70).  It  is  not  yet  appreciated  why  the  Mexican  coast  is 
passive  to  these  two  hurricanes. 

The  slow  oscillation  found  at  Key  West  and  St.  Petersburg  in  the 
Carla  simulation  is  also  present.  It  is  noticeable  again  at  all 
stations  from  Key  West  to  Appalachicola . 

The  maximum  computed  water  level  occurs  at  grid  point  (4,33) 
which  is  the  location  of  the  sample  point  for  the  Port  Isabel 
hydrograph.  This  simulated  peak  surge  is  of  the  order  of  one  meter 
higher  than  the  recorded  peak.  This  is  to  be  expected  considering 
the  difference  in  locations  of  the  simulated  and  actual  hydrographs 
(cf.  Fig.  29,  p.  72).  In  addition  to  reduction  of  the  magnitude  of 
wind  speed  at  the  actual  tide  station  by  South  Padre  Island,  the 
volume  responses  inside  the  Laguna  Madre  will  be  small  because  of  the 
constricted  opening.  These  factors  contribute  to  a  higher  computed 
water  level  at  this  station.  A  slower  retreat  in  the  observed 
hydrograph  is  due  possibly  to  the  trapping  of  water  inside  the 
embayment. 
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The  contours  of  modal  height  anomalies  and  surface  currents  at 
1200  GMT  9  August  are  shown  in  Figs.  47-49.  The  inertio-gravity 
waves  in  the  baroclinic  height  anomaly  (Fig.  48)  seem  to  be 
undetectable.  Examination  of  this  field  at  later  times  indicates 
that  only  one-half  wavelength  of  the  wake  oscillation  is  seen  in  Fig. 
48.  This  agrees  with  the  experimental  results  of  Chang  and  Anthes 
(1978)  showing  that  the  faster  the  storm  moves  the  longer  the 
wavelength  of  the  oscillation.  The  approximate  half  wavelength  of 
350  km  is  measured  in  this  hurricane  Allen  simulation. 

The  existence  of  f)G  in  the  simulated  hydrographs  for  Allen 
prompted  a  repeat  simulation  of  hurricane  Carla  with  the  added  code 
to  determine  tjq.  Figure  50  shows  the  resulting  time  sequence  of  jjg 
for  Carla.  The  amplitude  is  smaller  than  the  one  associated  with 
hurricane  Allen  and  the  averaged  period  is  about  24  h.  The  jjg  signal 
obtained  from  the  Carla  simulation  once  again  shows  up  simultaneouly 
in  the  individual  hydrographs  at  stations  around  the  Gulf. 

It  is  important  to  note  that  there  is  a  correlation  between  the 
»lG  signal  and  the  transport  through  both  Florida  and  Yucatan  Straits. 
Figures  51  and  52  show  the  time  series  of  the  volume  transport 
through  Florida  Strait  (FS),  Yucatan  Strait  (YS),  and  the  total 
differential  volume  transport  (unla^lled)  as  obtained  from  the 
hurricane  Carla  and  Allen  simulations.  The  striking  feature  of  the 
differential  transport  is  the  periodicity.  The  average  periods  of  24 
h  and  28  h  estimated  from  Figs.  51  and  52  are  exactly  the  same  as  the 
period  of  their  corresponding  ijg  signals.  In  both  cases, the  »iG 
signal  lags  the  net  transport  by  90*  in  phase.  The  first  maximum  of 
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Computed  barotropic  height  anomaly  field  (meters)  for  hurricane  Allen  at  1200  GMT  9  August 


Computed  surface  current  field  (cm/s)  for  hurricane  Allen  at  1200  GMT  9  August  1980 


5.0 


in  the  Gulf  of  Mexcio,  »jG,  obtained  from  the 
ea  level. 


Allen.  The  unlabelled  curve  is 


the  net  periodic  transport  produced  by  hurricane  Carla  is 
approximately  5  x  106  m^/s  as  compared  to  12.5  x  106  m3/s  generated 
by  hurricane  Allen.  The  t)q  signal  reaches  the  first  maximum  of  about 
0.1  m  for  hurricane  Carla  and  0.35  m  for  hurricane  Alien  which  are  in 
conformity  with  their  respective  net  periodic  transports.  These 
results  are  indicative  of  the  generation  of  by  the  net 
differential  port  forcing  (volume  transport). 
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CHAPTER  IV 


PARAMETRIC  STUDY 

The  purposes  of  the  parametric  study  were  to  obtain  responses  in 
the  Gulf  of  Mexico  to  different  forcing  from  hypothetical  storms  and 
to  investigate  those  cases  where  the  forerunner  surges  might  be 
generated. 

1.  Selection  of  paths 

A  total  of  5  paths,  designated  as  PATH1  to  PATH5,  for  which  the 
model  storms  would  traverse  the  Gulf  was  selected.  The  first  four 
paths  originated  in  the  Cayman  Sea  and  entered  the  Gulf  through 
Yucatan  Strait.  Sequentially,  the  locations  where  these  four  paths 
crossed  land  were  in  the  vicinity  of  Corpus  Christi  (PATH1),  Sabine 
Pass  ( PATH2 ) ,  Burrwood  ( PATH3 )  and  Apalachicola  (PATH4).  The  last 
path  ( PATH5 )  started  at  25°N  and  81.5°W  (overland  in  Florida)  and 
made  landfall  at  Corpus  Christi.  All  of  these  paths  were  great 
circles  as  shown  in  Fig.  S3.  The  first  and  second  track  (PATH1  and 
PATH2 )  were  similar  to  those  of  hurricanes  Carla  and  Allen.  The 
third  path  (PATH 3)  was  similar  to  hurricane  Camille’s  (1969)  track 
while  the  fourth  path  resembled  the  track  of  hurricane  Agnes  (1972). 
The  last  track  (PATHS)  which  is  rarely  observed  in  nature,  was 


selected  as  a  special  case. 


2.  Hypothetical  storms 


The  hypothetical  storms  were  assumed  to  be  characterized  by  the 
NWS23  parametric  model  with  time  invariant  radius  of  maximum  wind, 
translational  speed  and  central  pressure  deficit.  Thr^e  constant 
forward  speeds  of  15,  25  and  35  km/h  were  adopted.  The  slowest  one 
was  comparable  to  the  average  forward  speed  of  hurricane  Carla  while 
the  fastest  speed  was  approximately  the  translation  speed  of 
hurricane  Allen  when  it  was  in  the  Cayman  Sea. 

Two  radii  of  maximum  wind,  30  and  60  km,  were  selected  to 
account  for  the  scale  variation  of  the  atmospheric  forcing.  The 
radii  were  kept  constant  for  the  entire  simulation  period.  The  range 
of  radius  of  maximum  wind  for  Gulf  Coast  hurricanes  of  record  as 
summarized  in  Schwerdt  et  al.  (197 9)  was  approximately  10  to  60  km. 
This  range  of  radii  applies  to  those  storms  making  landfall  between 
Port  Isabel  and  Apalachicola,  which  covered  the  landfall  locations  of 
the  five  selected  paths.  The  large  radius  is  at  the  upper  limit  and 
the  other  is  about  the  average  of  the  observed  range.  The 
approximate  dimensions  of  a  grid  bloc*  are  25  km  x  27  km  and, 
therefore,  sets  a  lower  limit  on  the  radius  of  maximum  wind  which  can 
be  used.  Storms  with  a  radius  of  maximum  wind  smaller  than  the  grid 
size  would  not  be  resolved  and  the  associated  wind  field  would  be 
highly  distorted. 

The  last  characteristic  considered  was  the  pressure  drop  (AP 
=Pn-Pc) .  It  is  known,  however,  that  there  is  a  nearly  linear 
relationship  between  the  peak  surge  and  the  pressure  deficit,  other 
parameters  being  held  fixed  (Jelesnianski,  1972).  Therefore,  a 
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constant  value  of  80  mb  pressure  drop  was  assumed  in  most  of  these 
parametric  simulations.  Nevertheless,  one  of  the  sensitivity  tests, 
which  will  be  discussed  later,  was  designed  to  verify  this 
relationship. 

Table  1  summarizes  the  possible  combinations  of  these  parameters 
for  the  five  selected  paths  including  their  designated  hurricane 
names  for  further  discussions.  Runs  were  made  for  a  representative 
subset  of  these  model  hurricanes  (those  identified  by  an  *  in  the 
last  column) .  This  subset  allows  one  to  examine  responses  for 
different  paths  for  average  hurricane  parameters,  as  well  as  to 
determine  effects  of  forward  speed  and  effects  of  scale. 

3.  Simulation  procedure 

The  start  up  process  for  all  parametric  simulations  is  the  same 
as  in  the  simulations  of  hurricanes  of  record.  The  model 
calculations  start  at  0000  h  on  day  1  but  end  at  different  times 
depending  on  the  path  and  forward  speed.  However,  in  all  cases  the 
computations  proceeded  to  at  least  24  h  after  landfall. 

Wind  and  pressure  forcings  are  again  computed  at  each  time  step. 
Parameterizing  the  pressure  drop  instead  of  maximum  wind  speed 
(including  assuming  k=l)  simplifies  the  wind  stress  computations  by 
eliminating  linear  interpolations  in  space  (azimuth)  and  time  to 
obtain  k  and  Wffl.  Equation  (75)  is  used  to  compute  the  wind  speed 
which  is  then  augmented  by  the  forward  speed  to  account  for  the 
asymmetry  in  the  wind  fields.  The  only  information  required  to 
compute  the  forcing  is  the  position  of  the  storm  center  which  is 


Table  l.  Characteristics  and  designated  hurricane  names  for  the 
adopted  hypothetical  storms  for  parametric  study 


Radius  of  Forward 

Path  maximum  speed  Name 

wind  (km)  (km/h) 


* 


* 


30 


35 

25 

15 


HUR2 

HUR1 

HUR3 


given  at  6  h  intervals.  The  code  added  to  compute  r)G  is  retained. 


4.  Results  of  parametric  simulations 

The  simulated  responses  of  the  Gulf  to  the  synthetic  storms 
listed  in  Table  1  are  presented  in  this  section.  Discussions  of  the 
results  are  separated  into  three  parts.  First,  general  results 
common  to  all  simulations  are  briefly  discussed.  The  simulated 
hydrographs  obtained  from  each  simulation  and  other  results 
pertaining  to  the  development  of  forerunners  are  highlighted  in  the 
second  part  of  the  discussion.  The  last  part  contains  results  from 
related  simulations. 

a)  General  results 

The  inertio  -  gravity  waves  in  the  baroclinic  height  anomaly 
fields  are  found  in  all  parametric  runs.  Figures  54  through  57  show 
contours  of  baroclinic  height  anomalies  and  the  surface  current 
fields  obtained  from  the  HUR2  and  HUR3  simulations.  Increasing  the 
forward  speed  of  the  storm  resulted  in  increasing  the  wavelength  and 
decreasing  the  width  of  the  wake  as  seen  in  Figs.  54  and  56.  This  is 
in  agreement  with  the  results  of  Chang  and  Anthes  (1978)  and  Geisler 
(1970)  regarding  the  wavelength  and  the  width,  respectively,  of  the 
storm  generated  wake. 

A  dominant  long  period  oscillation  is  found  in  the  hydrographs 
from  all  stations  on  the  Florida  shelf.  Estimated  periods  of  this 
oscillation  are  on  the  order  of  4  -  6  days.  Marmorino  (1982) 
analyzed  sea  level  records  from  Key  West  to  Pensacola  recorded  during 
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Pig.  55.  Computed  surface  current  field  (cm/s)  at  0000  h  on  day  5  generated  by  model  storm  HUR2 


Computed  baroclinlc  height  anomaly  field  (meters)  at  0000  h  on  day  7  generated  by  model 


Computed  surface  current  field  (cra/s)  at  0000  h  on  day  7  generated  by  model  storm  HUR3 


January  -  April  1978.  Dominant  signals  were  found  in  the  6  day  band 
and  3  day  band.  However,  the  Gulf-wide  signal,  jjq,  can  also  be 
observed  in  the  simulated  hydrographs  from  stations  on  the  northern 
part  of  the  shelf.  Figures  58  and  55  show  hydrographs  from 
Apalachicola  and  Cedar  Key  obtained  from  the  HUR1  simulation.  It  is 
obvious  from  these  two  figures  that  other  modes  of  response  exist  on 
the  Florida  shelf. 

b)  Simulated  hydrographs 

The  following  discussions  are  ordered  according  to  the  selected 
paths  starting  from  PATH1  through  PATHS.  Within  each  path,  results 
obtained  from  the  storms  of  small  radius  are  presented  first  followed 
by  that  of  the  large  radius  storms . 

The  time  sequence  of  the  computed  water  levels  at  Galveston 
obtained  from  the  simulation  of  HUR1  is  shown  in  Fig.  60.  The 
initial  rise  of  water  level  which  reaches  a  maximum  of  0.22  m  at 
approximately  1300  h  on  day  4  is  well-defined.  Resurgence  after  the 
peak  surge  with  an  oscillation  period  of  about  6  h  is  also 
noticeable. 

The  initial  peak  of  water  level  in  the  Galveston  hydrograph  is 
in  phase  with  the  first  maximum  of  the  tjq  signal  as  shown  in  Fig.  61. 
The  first  extreme  of  nG,  which  is  also  the  maximum  peak,  occurs 
approximately  24  h  after  the  storm  entered  Yucatan  Strait.  The 
average  period  of  this  signal  is  30  h.  This  mode  is  not  only  present 
in  the  hydrographs  from  all  stations  but  also  accounts  for  almost  all 
of  the  response  in  the  southwest  sector  of  the  Gulf  as  shown  in  Figs. 


equence  of  water  levels  at  Apalachicola  obtained  from  the  simulation  of  HUR1 


Fig.  59.  Ti®2  sequence  of  water  levels  at  Cedar  Key  obtained 
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to.  Time  sequence  of  water  levels  at  Galveston  obtained  from  the  simulation  of  HUR1.  The  datum 
c;ean  sea  level . 


Hurricane 


Simulation  of  HUP. I .  The  solid  line  is  computed  by  averaging  water  levels  from  every  grid  point  in 
the  Gulf  at  each  time  step.  The  dashed  line  is  computed  from  the  continuity  equation.  The  datum  is 
r.can  sea  level. 


62  and  63.  The  dashed  lines  represent  rjG  and  the  solid  lines  are  the 
simulated  water  levels  from  the  indicated  station. 

The  maximum  surges  generated  by  HUR2  and  HUR3  imply  a  direct 
relationship  between  peak  surge  at  the  open  coast  and  forward  speed 
of  the  storm.  Jelesnianski  (1972)  proposed  a  correction  factor  for 
the  effect  of  storm  vector  motion  (track  direction  and  forward  speed 
at  landfall)  which  is  larger  for  a  faster  moving  storm  provided  that 
the  landfall  angles  are  the  same. 

Figures  64  and  65  show  the  simulated  hydrographs  at  Galveston 
from  th"i  HUR2  and  HUR3  simulations,  respectively.  The  initial  rise 
of  water  level  before  the  peak  surges  are  present  in  both  runs  which 
are  concurrent  with  the  first  maximum  of  their  corresponding  n q 
signals.  The  maximum  peaks  and  periods  of  riG  as  determined  from 
Figs.  66  and  67  are  0.22  m,  28  h  and  0.21  m  30  h,  respectively.  The 
time  lag,  &Qt  is  approximately  17  h  for  the  faster  storm  (HUR2)  and 
about  24  h  for  the  slower  storm  (HUR3). 

The  maximum  surges  produced  by  storms  of  large 
( HUR4 , HUR5 , HUR6 )  are  on  the  order  of  2  m  larger  than  those 
corresponding  to  the  small  storm  (RUR1,HUR2,HUR3)  simulations. 

Except  for  a  larger  percentage  of  increase  in  the  peak  surge  at  the 
open  coast,  the  results  are  in  qualitative  agreement  with 
Jelesnianski  (1972)  in  which  a  very  simple  bathymetry  was  employed. 

The  hydrograph  at  Galveston  from  the  KUR4  simulation  is  shown  in 
Fig.  68.  The  time  of  the  peak  surge  and  the  time  at  which  the 
initial  rise  of  water  level  reaches  the  maximum  are  exactly  the  same 
as  in  the  KUR1  simulation.  The  resurgence  oscillation  is  concealed 


Pig.  62.  Time  sequence?  of  water  levels  at  Port  Isabel  (solid)  and  n.Q  (dashed)  obtained  from  the 
simulation  of  HUR1.  The  datura  is  mean  sea  level. 


ig.  63.  Time  sequences  of  water  levels  at  Proyreso  (solid)  and  rjG  (dashed)  obtained  from  the 
simulation  of  HUR1.  The  datum  is  mean  sea  level. 
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65.  Time  sequence  of  water  levels  at  Galveston  obtained  from  the  simulation  of  HUR3.  The  datum 
lean  sea  level. 


Fig.  67.  Time  sequences  of  the  average  water  levels  in  the  Gulf  of  Mexico,  r? q,  obtained  from  the 
simulation  of  HUR3 .  The  solid  line  is  computed  by  averaging  water  levels  from  every  grid  point  in 
the  Gulf  at  each  time  step.  The  dashed  line  is  computed  from  the  continuity  equation.  The  datum  is 
mean  sea  level. 


sequence  of  water  levels  at  Galveston  obtained  from  the  simulation  of  HUR4.  The  datum 


by  a  high  peak  surge  but  still  observable.  Figure  69  shows  the  time 
sequence  of  the  Gulf-wide  oscillation,  tjg,  generated  by  HUR4.  The 
average  period  and  6G  are  30  h  and  24  h,  respectively.  The  maximum 
peak  is  0.42  m,  about  a  factor  of  2  larger  than  in  the  HUR1 
simulation.  It  is  interesting  to  note  that  the  largest  difference  in 
the  extremes  of  rjG  from  HUR1  and  HUR4  occur  at  their  first  maximum. 
During  later  stages  the  differences  reduce  to  a  few  centimeters.  The 
presence  of  jjg  in  the  simulated  hydrographs  from  the  other  stations 
is  preserved.  However,  this  oscillation  is  no  longer  coincident  with 
the  water  levels  at  the  individual  stations  as  in  the  previous  cases. 
Figures  70  and  71  clearly  demonstrate  the  departure  between  the  two 
curves  at  Port  Isabel  and  Progreso. 

A  peak  surge  of  6.5  m  at  the  coast  is  generated  by  a  fast  moving 
storm  of  large  Rmav(HUR5) .  The  slow  moving  storm  of  the  same  size 
( HUR6 )  however,  produces  a  smaller  4.71  m  peak  surge.  Increasing  the 
maximum  surge  at  the  coast  with  increasing  forward  speed  is  the  same 
as  that  obtained  from  simulations  of  small  storms  as  discussed  above. 
Pigures  72  and  73  display  the  computed  hydrographs  at  Galveston 
obtained  from  the  HUR5  and  HUR6  simulations,  respectively.  Both 
hydrographs  have  an  initial  rise  of  water  level  prior  to  the  peak 
surge  that  again  matches  the  first  peak  of  their  corresponding  ijq 
series.  The  time  sequence  of  the  nG  signal  for  KUR5  and  HUR6  are 

shown  in  Figs.  74  and  75.  The  maximum  peaks  of  ijg  (0.4  m  for  HUR5 

and  0.38  m  for  HUR6)  vary  slightly  with  the  forward  speed  hut  both 
are  approximately  twice  as  large  compared  with  those  produced  by  the 

small  storms  with  the  same  forward  speed,  i.e.,  HUR2  and  HUR3, 
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simulation  of  HUR4.  The  solid  line  is  computed  by  averaging  water  levels  from  every  grid  point  in 
the  Gulf  at  each  time  step.  The  dashed  line  is  computed  from  the  continuity  equation.  The  datum  is 
mean  sea  level. 
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Fig.  70.  Time  sequences  of  water  levels  at  Port  Isabel  (solid)  and  rjG  (dashed)  obtained  from  the 
simulation  of  HUR4.  The  datura  is  mean  sea  level. 


ig .  71.  Time  sequences  of  water  levels  at  Pregre 
simulation  of  IIUR4  .  The  datum  is  mean  sea  level. 


72.  Time  sequence  of  water  levels  at  Galveston  obtained  from  the  simulation  of  HIJR5.  The  datum 
mean  sea  level. 


Ig.  73.  Time  sequence  of  water  levels  at  Galveston  obtained  from  the  simulation  of  HUR6.  The  datum 
is  mean  sea  level. 
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lg,  74.  Time  sequences  of  the  average  water  levels  In  the  Gulf  of  Mexico,  tjq,  obtained  from  the 
simulation  of  HUR5,  The  solid  line  is  computed  by  averaging  water  levels  from  every  grid  point  in 
the  Gulf  at  each  time  step.  The  dashed  line  is  computed  from  the  continuity  equation.  The  datum  is 
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respectively.  Comparisons  of  the  tjg  hydrographs  between  HUR5  and 
HUR2  and  also  between  HUR6  and  HUR3,  reveal  that  the  largest 
differences  occur  in  the  early  stages  of  the  simulations.  This  is 
the  same  result  obtained  from  similar  comparisons  between  the  HUR1 
and  HUR4  Gulf-wide  modes. 

Important  results  from  the  PATH1  simulations  are  summarized  as 
follows.  The  initial  rise  of  water  level  before  the  peak  surge  at 
Galveston  exists  in  all  cases.  This  rise  matches  the  first  peak  of 
the  Gulf-wide  oscillation,  rjG,  which  occurs  simultaneously  around  the 
Gulf.  Periods  of  r)G  vary  between  28  h  to  30  h.  The  effects  of  storm 
size  (Rjnax)  are  to  increase  both  the  peak  surge  at  the  coast  and  the 
magnitude  of  tjq-  On  the  contrary,  the  effect  of  the  storm  forward 
speed  is  to  increase  only  the  peak  surge  at  the  coast  (for  the  same 
size  storms).  Accordingly,  for  the  remaining  storm  tracks  ( PATH2 
through  PATH5)  only  one  simulation  for  each  path  using  a  storm  of  30 
km  Rjj^  and  25  km/h  forward  speed  was  made. 

A  peak  surge  of  4.4  m  at  grid  point  (23,46)  near  5W  Pass, 
Louisiana,  is  obtained  from  the  PATH2  ( HUR7 )  simulation.  The  nearest 
station  to  this  point  is  Grand  Isle  which  is  situated  to  the  east. 

The  time  history  of  the  computed  water  levels  at  this  station  as 
shown  in  Fig.  76.  The  highest  computed  water  level  at  Grand  Isle 
occurs  a  few  hours  earlier  than  the  true  peak  surge  (at  grid  point 
23c46)  while  the  first  maximum  in  this  hydrograph  coincides  with  the 
first  extreme  of  »jg.  The  sequential  occurrence  of  the  two  peaks 
resemble  the  situation  at  Galveston  in  the  PATH1  simulations.  The  ijg 
hytfrograph  shown  in  Fig.  77  indicates  the  first  maximum  of  0.2  m  with 
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Fig.  76.  Time  sequence  of  water  levels  at  Grand  Isle  obtained  from  the  simulation  of  HUR7 .  The  datum 
is  mean  sea  level. 
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Fig.  77.  Tima  sequences  of  the  average  water  levels  in  the  Gulf  of  Mexico,  t)q,  obtained  from  the 
simulation  of  HUR7 .  The  solid  line  is  computed  by  averaging  water  levels  from  every  grid  point  in 
the  Gulf  at  each  time  step.  The  dashed  line  is  computed  from  the  continuity  equation.  The  datum  is 
mean  sea  level. 


5G  of  approximately  19  h.  The  estimated  period  is  28  h.  This  Gulf- 
wide  oscillation  resumes  its  dominant  role  in  the  computed 
hydrographs  on  the  left  side  of  the  track  from  Port  Isabel  to 
Progreso. 

The  model  storm  HUR13  which  traversed  PATH3  produced  a  peak 
surge  of  5.81  m  near  Burrwood  (grid  point  37,48).  The  hydrograph 
from  Pensacola  shown  in  Fig.  78  indicates  a  maximum  water  level  of 
less  than  one  meter.  The  difference  of  almost  5  m  in  maximum  surge 
between  Pensacola  and  Burrwood  is  probably  due  to  the  very  narrow 
shelf  width  at  Pensacola.  The  initial  rise  of  water  levels  that 
match  the  first  maximum  of  tjg  is  noticeable.  The  period,  maximum 
peak  and  time  lag,  SG,  determined  from  Fig.  79  are  26  h,  0.2  m  and  11 
h,  respectively.  The  t iq  signal  remains  detectable  in  the  hydrographs 
from  stations  around  the  Gulf  and  again  is  the  major  part  of  the 
response  at  stations  on  the  west  and  southwest  coast. 

The  storm  moving  due  north  along  PATH4  at  25  km/h  with  a  30  km 
radius  of  maximum  wind  (HUR19)  generated  a  maximum  surge  of  6.22  m 
one  grid  block  to  the  east  of  the  Apalachicola  sampling  point. 
However,  the  simulated  hydrographs  show  higher  water  levels  at  Cedar 
Key  than  at  Apalachicola.  The  sampling  point  for  Apalachicola  is 
only  half  a  grid  block  away  from  the  path.  Consequently,  this 
station  lies  inside  the  radial  distance  between  the  storm  center  and 
the  maximum  wind  location  and  is  hence  subjected  to  a  smaller  wind 
stress.  The  time  sequence  of  the  computed  water  level  at  Cedar  Key 
is  shown  in  Fig.  80.  The  initial  rise  of  water  level  prior  to  the 
peak  surge  is  almost  undetectable.  The  tjg  hydrograph,  Fig.  81, 
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e  sequence  of  water  levels  at  Pensacola  obtained  from  the  simulation  of  HUR13 .  The  datum 
level . 


urn 


of  HUR13 .  The  solid  line  is  computed  by  averaging  water  levels  from  every  grid  point  in 
t  each  time  step.  The  dashed  line  is  computed  from  the  continuity  equation.  The  datum  is 
evel . 


indicates  the  first  maximum  of  only  a  few  centimeters.  The  estimated 
period  of  jjg  and  SG  are  27  h  and  12  h,  respectively.  In  spite  of  its 
small  amplitude,  r;G  still  dominates  the  response  in  the  hydrographs 
from  Pensacola  counterclockwise  to  Progreso.  The  striking  result 
however,  is  the  generation  of  a  12  h  period  oscillation  after  the 
peak  surge.  This  mode  is  also  observed  in  the  Apalachicola  and  St. 
Petersburg  hydrographs. 

The  very  small  initial  rise  of  water  levels  along  the  Florida 
shelf  prior  to  the  peak  surge,  together  with  the  generation  of  the  12 
h  period  oscillation  motivated  three  additional  storm  runs  for  PATH 4 . 
All  have  large  but  take  three  different  forward  speeds. 

The  medium  speed  and  large  radius  storm  (HUR22)  produced  a  8.53 
m  maximum  surge  to  the  north  of  Cedar  Key.  The  simulated  hydrograph 
at  Cedar  Key,  Fig.  82,  shows  a  peak  surge  of  more  than  5  m  followed 
by  a  strong  12  h  period  oscillation.  The  presence  of  the  initial 
rise  of  water  level  is  hardly  seen.  Figure  83  reveals  that  a  first 
maximum  of  less  than  0.1  m  occurs  at  approximately  180?  h  on  day  4. 

At  this  time  the  water  level  at  Cedar  Key  is  decreasing  and  reaches  a 
minimum  at  1200  h  on  day  4.  The  stipulated  path  and  a  large  radius 
of  maximum  wind  resulted  in  larger  draw  down  of  water  level  in  the 
early  stages  due  to  stronger  offshore  directed  wind  preceding  the 
storm  as  compared  to  the  simulation  of  HUR19.  The  hydrograph  of  tjG 
shows  a  maximum  peak  of  0.24  m  at  2300  h  on  day  4.  An  estimated 
period  of  29  h  is  determined  from  Fig.  83.  Once  again,  the 
hydrographs  from  stations  on  the  left  side  of  the  track  from 
Pensacola  to  Progreso  seem  to  depart  from  nG  signal,  the  situation 
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Time  sequence  of  water  levels  at  Cedar  Key  obtained  from  the  simulation 
tsean  sea  level. 


Fig.  CJ.  Time  sequences  of  the  average  water  levels  in  the  Gulf  of  Mexico,  r?G,  obtained  from  the 
simulation  of  BUR22,  The  solid  line  is  computed  by  averaging  water  levels  from  every  grid  point  in 
the  Gulf  at  each  time  step.  The  dashed  line  is  cofiiputed  from  the  continuity  equation.  The  datum  is 


found  in  simulating  large  storms  on  PATH1. 

A  maximum  surge  of  8.52  m  occurring  at  the  same  location  as  in 
the  HUR22  run  was  generated  by  HUR23 .  The  time  sequence  of  water 
levels  at  Cedar  Key  is  shown  in  Fig.  84.  Due  to  a  fast  speed  and  a 
relatively  short  distance  across  the  Gulf,  the  maximum  surge  at  Cedar 
Key  occurs  so  early  that  it  is  not  possible  to  identify  the  existence 
of  an  initial  rise  of  water  levels.  The  12  h  mode  however,  remains 
noticeable.  The  r )G  signal  (Fig.  85)  reaches  the  first  peak  of  0.2  m 
at  1400  h  on  day  3,  which  is  approximately  9  h  after  the  storm  passed 
Yucatan  Strait.  The  highest  peak  of  rjG  is  0.23  m  and  the  period  is 
estimated  at  25  h. 

Reducing  the  forward  speed  of  the  storm  to  15  km/h  (HUR24), 
decreases  the  maximum  surge  to  7.32  m.  Figure  86  illustrates  two 
negligible  water  level  maxima  before  the  peak  surge  at  Cedar  Key. 

The  12  h  mode  after  the  peak  surge  is  again  excited.  The  26  h  period 
of  the  small  amplitude  rjG  signal  is  shown  in  Pig.  87.  It  is 
interesting  but  not  yet  understood  why  the  first  maximum  of  ijg 
occurred  before  the  storm  center  entered  the  Gulf.  The  6q  determined 
from  the  second  maximum  of  gG  is  approximately  20  h. 

These  four  runs  on  PATH4  repeatedly  gave  the  same  results 
regarding  the  effects  of  forward  speed  and  radius  of  maximum  wind. 

The  notable  result  is  the  excitation  of  the  12  h  mode  on  the  Florida 
shelf.  It  is  important  to  emphasise  that  this  model  admitted  no 
tidal  forcing.  Reid  and  Whitaker  (1981)  found  a  near  resonant 
response  on  the  Florida  shelf,  with  the  greatest  signal  near  Cedar 
Key,  to  direct  forcing  by  the  M,  tide  potential.  Figure  88  (Ictuye 


fig.  84.  Time  sequence  of  Mater  levels  at  Cedar  Key  obtained  from  the  simulation  of  HUR23 .  The  datum 


66.  T'rae  sequence  of  water  levels  at  Cedar  Key  obtained  from  the  simulation 
mean  sea  level. 


et  al.,  1973)  shows  the  recorded  water  level  with  tidal  signal 


removed  at  Cedar  Key  during  the  passage  of  hurricane  Agnes  in  1972. 
The  12  h  period  oscillation  is  obvious. 

The  maximum  surges  on  the  order  of  7-8  m  on  the  Florida  shelf 
found  in  the  PATH 4  simulations  are  probably  too  large.  Reid  and 
Whitaker  (1981)  showed  that  the  response  on  the  Florida  shelf  to 
direct  forcing  by  the  M2  tide  potential  was  very  sensitive  to  the 
local  friction  coefficient.  Therefore,  the  simulation  of  HUR23  was 
repeated  with  a  friction  coefficient  of  7.5  x  10-4  m/s.  This 
coefficient  is  applied  to  the  entire  computing  domain  since  only  the 
response  on  the  Florida  shelf  is  the  primary  concern  in  this 
simulation.  The  maximum  surge  generated  (6.23  m)  is  decreased  by 
about  25  %  from  that  generated  by  the  same  storm  simulation  with  the 
smaller  friction  coefficient.  The  hydrographs  from  stations  on  the 
Florida  shelf  (not  shown)  clearly  demonstrated  the  effect  of  larger 
friction.  Examination  of  the  remaining  hydrographs  reveals  no 
significant  changes  in  water  levels  due  to  increasing  friction,  even 
on  the  Texas-Louisiana  shelf  region.  Presumably  these  stations  are 
located  too  far  to  the  left  of  the  storm  track  such  that  the 
corresponding  currents  and  bottom  stress  are  small. 

The  PATHS  storm,  HUR25,  produced  a  5.4  m  maximum  surge  near 
Oalveston,  Fig.  89.  The  sharp  drop  of  water  levels  before  the  peak 
surge  is  caused  by  a  direct  offshore  wind  stress.  The  maximum  peak 
of  nG  (Fig.  90)  is  less  than  0.1  o  .The  estimated  period  of  »|q  for 


this  run  is  32  h. 


me  sequence  of  water  levels  at  Galveston  obtained  from  the  simulation  of  HUR25 .  The  datum 
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of  HOR25.  The  solid  line  is  computed  by  averaging  water  levels  from  every  grid  point  in 
each  time  step.  The  dashed  line  is  computed  from  the  continuity  equation.  The  datum  is 


The  relationship  between  the  magnitude  of  j)(j  and  the  net  volume 
volume  transport  through  the  ports  discussed  in  the  simulation  of 
hurricane  Allen  is  reinvestigated.  Figures  91  through  93  show  the 
barotropic  volume  transport  through  the  Florida  port  and  Yucatan 
Strait,  and  the  difference  of  the  two  components  obtained  from  the 
HUR1,  HUR4  and  HUR25  simulations,  respectively.  It  can  be  seen  from 
the  figures  that  the  larger  the  net  volume  transport  through  the 
ports  the  larger  the  magnitude  of  r>G.  This  result  should  be  expected 
based  on  continuity.  The  first  maximum  of  the  total  volume  transport 
for  KUR4  is  approximately  10  x  106  m3/s  which  is  about  a  factor  of  2 
larger  than  that  produced  by  HUR1.  The  same  ratio  of  maximum  i?G 
obtained  from  the  two  storms  provide  quantitative  evidence  for 
relating  tjg  to  the  net  periodic  transport.  Note  that  the  volume 
transport  from  the  HUR25  simulation  and  the  associated  gG  is  very 
small.  The  initial  position  of  this  storm  and  associated  wind  field 
gives  smaller  flows  through  the  two  openings  as  compared  to  HUR1  and 
HUR4 . 

The  average  period  of  the  net  periodic  volume  transport  in  each 
case  is  about  the  same  as  the  period  of  the  corresponding  rjG  signals. 

O 

The  nG  signal  however,  is  approximately  90  out  of  phase  (lag)  with 
the  total  periodic  transport.  This  result  resembles  the  uninodal 
seich  in  a  one  dimensional  channel.  However,  the  ubiquity  of  the  i|q 
signal  in  the  hydrographs  around  the  Gulf  implies  that  the  net 
periodic  transport  acts  like  a  single  port  forcing  and  excites  the 
entire  Gulf  proper  to  oscillate  with  a  node  at  the  port.  In  other 
words,  the  nG  oscillation  is  a  Helmholtz  mode  excited  by  the  net 
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sequences  of  volume  transport  through  Florida  (FS)  and  Yucatan  (YS)  Straits  obtained 
ation  of  HUR1.  The  unlabelled  curve  is  the  net  differential  volume  transport. 


Into  sequences  c£  volume  transport  through  Florida  (FS)  and  Yucatan  (YS)  Straits  obtained 
insulation  of  HUR4 .  The  unlabelled  curve  is  the  net  differential  volume  transport. 
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sequences  of  volume  transport  through  Florida  (FS)  and  Yucatan  (YS)  Straits  obtained 
at  ion  of  HUR2‘>.  The  unlabelled  curve  is  the  net  differential  volume  transport. 


periodic  transport  through  the  ports. 

Table  2  summarizes  the  peak  surge  at  the  digitized  coastline  for 
each  storm.  Locations  where  the  peaks  occur  are  given  in  terms  of 
the  grid  point  (I,J).  The  peaks  and  periods  (TQ)  of  the  tjg  signal 
obtained  from  each  case  are  given  in  columns  7,8  and  13,14.  The  time 
lag  of  the  first  maximum  of  jjg  after  the  center  of  the  storms  entered 
Yucatan  Strait,  6q,  is  also  given  in  Table  2. 

5.  Long  period  variation  of  jjg 

a)  The  3.4  day  volume  mode 

The  j]q  signals  obtained  from  all  the  synthetic  storm  simulations 
have  a  long  period  variation  superimposed  on  their  average  28  h 
period.  This  slow  variation  is  not  readily  observable  except  for  the 
one  generated  by  HUR5 .  Therefore,  model  HUR5  was  used  to  investigate 
this  long  period  variation  of  the  Helmholtz  mode  (>1q).  This 
simulation  used  the  same  storm  characteristics  of  the  original  HUR5 
hurricane.  However,  forcing  was  allowed  only  over  the  first  3  days 
22  h  of  the  simulation.  The  calculations  were  continued  for  a  total 
of  21  days  to  provide  a  clear  history  of  the  free  modes  in  the  Gulf 
and  Cayman  Sea.  Time  histories  of  the  mean  external  height  anomalies 
of  the  two  basins  over  the  full  simulation  period  are  shown  in  Fig. 
94.  The  arrow  indicates  when  all  forcing  was  set  to  zero.  At  this 
time  HUR5  was  centered  at  25.44*N  and  90.66*W. 

The  most  striking  feature  of  this  simulation  is  the  very  long 
code  which  dominates  the  Cayman  response.  The  period  of  this  mode  is 
approximately  3.4  days.  The  half -cycle  means  of  Qq  over  the  last  7 
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days  of  tha  simulation  correspond  very  closely  to  the  3.4  day  mode  of 
the  Cayman  Sea.  This  implies  a  complicated  response  of  the  Gulf- 
Cayman  Sea  system  which  oscillates  in  unison  at  a  period  of  3.4  days. 
Simultaneously  the  Gulf  of  Mexico  exhibits  a  28  h  volume  mode 
superimposed  on  the  long  period  oscillation.  The  latter  mode  can 
clearly  affect  the  time-dependent  amplitudes  assigned  to  riq- 

The  intriging  questions  of  course  concern  the  nature  of  the  3.4 
day  oscillation.  Definitive  answers  are  outside  the  scope  of  this 
research.  However,  the  influence  of  the  earth's  rotation  could  be  so 
readily  evaluated  that  it  was  decided  to  repeat  the  above  simulation 
with  the  Coriolis  acceleration  set  to  zero.  These  results  are  shown 
in  Fig.  95. 

The  3.4  day  mode  was  not  excited  in  the  non-rotating  Gulf  and 
Cayman  basins .  Inspection  of  the  -accompanying  digital  means  did  not 
reveal  a  discernible  3  day  mode.  Moreover,  tne  digital  time  series 
show  that  the  28  h  volume  mode  in  the  Gulf  was  reduced  to  amplitudes 
on  the  order  of  a  centimeter. 

Notice  that  the  paired  histories  of  water  level  with  (cf.  Fig. 
94,  p.  158)  and  without  (Fig.  95)  rotation  are  different  over  the 
entire  period  of  forcing  (first  3  days  22  h).  This  suggests  that  the 
excitation  mechanism  for  the  Helmholtz  mode  in  the  Gulf  and  the  3.4 
day  Gulf-Cayman  system  r.ode  is  contingent  on  the  earth's  rotation. 


Hydrograph  (ETAG) 


Fig.  95.  Time  sequences  of  the  averaged  water  level  in  the  Gulf  of  Mexico  (solid)  and  in  the  Cayman 
Sea  (dashed)  obtained  from  the  long  term  simulation  of  HUR5  with  no  rotation.  The  datum  is  mean  sea 
level . 


b)  The  6.5  day  tilt  mode 


Time  histories  of  volume  transport  obtained  from  the  21  day 
simulation  of  HUR5  with  rotation  are  presented  in  Fig.  96.  A 
prominent  feature  is  the  6.5  day  period  oscillation  for  the 
transports  through  both  Florida  and  Yucatan  Straits. 

The  water  level  from  this  long  simulation  was  sampled  at  12 
stations,  six  of  which  are  the  same  as  those  given  in  Fig.  12.  The 
remaining  six  stations  are  Panama  City,  Tampico,  Coatzacoalcas,  Dimas 
(northern  coast  of  Cuba, ,  Central  Gulf,  Eastern  Gulf  and  Western 
Gulf.  The  hydrograph  at  Key  West  (Fig.  97)  shows  a  strong  6.5  day 
mode  with  an  average  amplitude  of  0.50  m.  Superimposed  on  this  long 
period  mode  is  the  28  h  volume  mode.  Careful  inspection  of  the 
remaining  hydrographs  indicates  the  presence  of  the  6.5  day 
oscillation  but  with  a  much  smaller  amplitude  as  compared  to  that  in 
the  Key  West  hydrograph.  Hydrographs  at  Galveston  (Fig.  98)  and 
Dimas  (Fig.  99)  show  that  the  water  levels  at  the  two  stations  are 
180°  out-of -phase.  For  example,  at  1600  h  on  day  8  the  water  level 
at  Galveston  is  0.05  m  while  at  Dimas  it  is  -0.05  m.  The  three  deep¬ 
water  hydrographs  are  near  zero  at  this  time  (note  that  the  Central 
Gulf  station  is  approximately  midway  between  Dimas  and  Galveston) . 
This  0.1  m  differ  itiax  in  water  surface  elevation  is  possibly  a 
geostrophic  tilt  since  all  the  forcings  were  set  to  zero  well  before 
this  time.  Large  space  and  time  scales  and  a  small  surface 
expression  are  indicative  of  a  vorticity  (or  quasi-geostrophic  tilt) 
mode. 
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g.  97.  Time  sequence  of  w 
he  datum  is  mean  sea  leve? 


Hurricane  5 

Dimas  (Cuba)  Hydrograph 
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HUR5.  The  datum  is  mean  sea  level 


A  simple  geostrophic  calculation  using  an  instantaneous  volume 
transport  through  F?ftida  Strait  at  1600  h  on  day  8  of  30  x  106  m^/s 
and  an  average  depth  across  the  port  of  950  m  yield  a  northward 
surface  gradient  from  Cuba  of  0.2  m.  A  similar  computation  for  the 
Yucatan  Strait  yields  a  surface  gradient  toward  Campeche  Bank  of  only 
0.07  m.  The  surface  height  field  is  not  available  for  comparison. 
However,  using  the  mean  water  levels  at  Dimas  and  Key  West  as  a  rough 
representation  of  the  north-south  surface  component  of  gradient 
reveals  that  the  tilting  of  water  surface  across  Florida  Strait  is  on 
the  order  of  0.5  m.  This  result  is  much  larger  than  that  obtained 
from  the  simplified  geostrophic  calculation.  It  is  quite  possible 
that  the  current  speed  is  not  uniform  across  the  ports.  An  estimate 
of  the  effect  of  current  shear  on  the  cross-port  surface  gradient  by 
numerically  integrating  the  geostrophic  equation  across  the  port  was 
made.  This  integration  is  given  by  the  relation 

=  Tj-j  ♦  (Ayfu/gDj^  ,  (79: 

where  U  is  the  transport  per  unit  width  for  each  grid  block  across 
the  Florida  Strait  and  Dj^  is  the  average  depth  for  each  block. 

Upon  assuming  uniform  distribution  of  transport  across  the  port,  U  is 
Q/W  where  Q  is  the  total  volume  transport  and  W  is  the  port  width. 
Employing  the  same  value  of  Q  as  in  the  previous  calculation  and 
taking  the  water  level  at  Cuba  (south  of  Key  West)  as  zero  yields  a 
surface  gradient  between  Key  West  and  Cuba  of  0.5  o.  For  the  Yucatan 
Strait  the  tilting  of  water  surface  between  Cuba  and  the  east  coast 
of  Yucatan  obtained  from  this  computation  is  0.26  m.  Even  though  Q 


is  the  same  at  both  ports,  Yucatan  channel  is  much  deeper  than  the 
Florida  port  so  that  the  computed  surface  gradient  is  smaller.  There 
is  no  hydrograph  on  the  east  coast  of  Yucatan  for  comparison  in  this 
case. 

6.  Results  from  related  simulations 

a)  Variation  of  pressure  drop 

All  the  synthetic  storms  employed  a  constant  pressure  drop  of  80 
mb.  As  a  sensitivity  test  of  the  model,  simulations  of  HUR5  with  a 
40  mb  ( HUR5W )  and  a  120  mb  (HUR5S)  pressure  drop  were  made.  Maximum 
surges  of  3.59  m  and  9.35  m  were  generated  in  the  simulations  of 
HUR5W  and  HUR5S,  respectively.  Figure  100  displays  the  maximum  surge 
as  a  function  of  pressure  drop  obtained  from  all  simulations  of  HUR5 . 
This  result  is  consistent  with  Jelesnianski ' s  (1972)  inference  that 
the  peak  surge  is  almost  a  linear  function  of  the  pressure  drop. 

Plots  of  obtained  from  the  HUR5W  and  HUR5S  simulations  are 
shown  in  Figs.  101  and  102,  respectively.  The  maximum  peak  and  the 
period  of  this  oscillation  are  0.2  m  and  27  h  for  the  weak  storm  and 
0.6  m  and  27  h  for  the  intense  storm. 

Note  that  the  HUR5W  simulation  yields  a  larger  magnitude  of  i )G 
compared  with  HUR25  despite  its  smaller  pressure  drop.  This  implies 
that  a  major  factor  governing  the  magnitude  of  the  from  these 
simulations  is  the  differential  wind-driven  volume  transport  through 
Florida  port  and  Yucatan  Strait. 


Pressure  Drop  (mb) 


Fig.  100.  Relationship  between  the  maximum  peak  surge  at  the  coast 
and  the  central  pressure  deficit  obtained  from  three  simulations  of 
HUR5  with  40  mb,  80  mb  and  120  mb  pressure  drop. 
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continuity  equation.  The  datum  is  mean  sea  level. 


lg.  iuz.  Time  sequences  of  the  average  water  levels  in  the  Gulf  of  Mexico,  rjG,  obtained  from  the 
simulation  of  HUR5  with  120  mb  pressure  drop  (HUR5S) .  The  solid  line  is  computed  by  averaging  water 
levels  from  every  grid  point  in  the  Gulf  at  each  time  step.  The  dashed  line  is  computed  from  the 
continuity  equation.  The  datum  is  mean  sea  level. 


b)  Barotropic  model 


It  is  important  to  ascertain  the  effects  of  the  baroclinic 
response  in  deep  water  on  the  maximum  surge  at  the  coast  since  both 
barotropic  and  baroclinic  modes  are  coupled  through  the  gradient 
terms,  VCH^/D).  Therefore,  an  additional  simulation  of  HUR5  was 
repeated  which  excluded  the  baroclinic  mode. 

A  maximum  peak  surge  of  6.42  m  at  grid  point  (14,46)  was 
obtained  from  this  simulation.  This  peak  is  only  eight  tenths  of  a 
centimeter  lower  than  that  obtained  from  the  two  mode  simulation  of 
this  storm.  Comparisons  of  the  simulated  hydrographs  from  the 
barotropic  model  and  the  original  HUR5  reveal  that  the  baroclinic 
response  does  not  produce  a  significant  contribution  to  the 
barotropic  response  on  the  shelf  except  at  Key  West,  Naples  and  St. 
Petersburg.  The  solid  and  dashed  lines  in  Pig.  103  and  104, 
respectively,  represent  the  computed  hydrographs  obtained  from  the 
two  mode  and  pure  barotropic  mode  simulations  at  the  indicated 
stations.  The  reason  that  this  effect  is  visible  on  the  Florida 
shelf  rather  than  other  locations  may  be  due  to  a  larger  baroclinic 
transport  (through  the  Florida  strait)  and  a  stronger  coupling  as  a 
result  of  steeper  slopes  at  the  Florida  shelf  break. 

As  expected,  the  Gulf-wide  oscillation  as  shown  in  Fig.  105 
obtained  from  this  simulation  is  only  slightly  different  from  that 
obtained  from  the  two  mode  run. 


rinse  sequences  of  water  levels  at  Cedar  Key  obtained  from  the  simulations  of  HUR5  (dashed) 
-e  barotrcpic  model  of  HUR5  (solid).  The  datum  is  mean  sea  level. 


Petersburg  Hydrograph 


ig.  104.  Time  sequences  of  water  levels  at  St.  Petersburg  obtained  from  the  simulations  of  HUR5 
(dashed)  and  the  pure  barotropic  model  of  HUR5  (solid).  The  datum  is  mean  sea  level. 
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ig •  1Q5.  Time  sequences  of  the  average  water  levels  in  the  Gulf  of  Mexico,  r?G,  obtained  from  the 
simulation  of  the  pure  barotropic  model  of  HUR5 .  The  solid  line  is  computed  by  averaging  water 
levels  from  every  grid  point  in  the  Gulf  at  each  time  step.  The  dashed  line  is  computed  from  the 
continuity  equation.  The  datum  is  mean  sea  level. 


c)  Radiation  boundary  condition 


The  effects  of  the  open  boundary  condition  employed  in  this 
study  on  rjG  was  of  primary  concern.  In  order  to  evaluate  this 
effect,  a  pure  radiation  boundary  condition  was  imposed  in  the 
simulations  of  HUR5,  HUR23  and  hurricane  Allen.  The  names  of  these 
storms  modified  by  (R)  are  used  in  the  following  discussions  to 
distinguish  them  from  their  corresponding  original  simulations. 

Figures  106  through  108  show  the  time  histories  of  tjg  obtained 
from  the  simulations  of  HUR5(R),  HUR23(R)  and  hurricane  Allen(R), 
respectively.  Only  the  r>Q  signal  obtained  by  averaging  water  levels 
in  the  Gulf  is  shown  in  the  figures  (note  the  changes  in  height 
scale) .  The  radiation  condition  at  the  open  boundaries  effectively 
radiates  tjg  out  of  the  Gulf.  The  peaks  of  tjg,  sequentially,  are  0.28 
m,  0.14  m  and  0.22  m  for  HUR5(R),  HUR23(R)  and  hurricane  Allen(R). 

The  estimated  e-folding  times  of  the  damping  rate  are  34  h  for 
HUR5(R) ,  24  h  for  HUR23(R)  and  less  than  10  h  for  hurricane  Allen (R) . 
The  radiation  condition  not  only  drastically  damps  the  rjG  signal  but 
also  changes  their  periods.  Reid  and  Whitaker  (1981)  experimentally 
determined  the  damping  rate  of  the  volume  mode  in  the  GOMT  model  to 
be  2.81  x  10”®  s-1  or  an  e-folding  time  of  2.6  days.  The  radiation 
condition  employed  in  the  GOMT  model  has  a  complex  admittance 
coefficient  in  which  the  imaginary  part  governed  the  effective  added 
mass  of  adjoining  seas.  The  radiation  condition  employed  in  this 
study  is  a  special  case  where  only  the  real  part  of  the  admittance  is 
considered.  Exclusion  of  the  imaginary  part  results  in  a  much  larger 
damping  rate  of  the  volume  mode  as  found  in  this  study.  However,  the 
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simulation  of  HUR5  with  radiation  condition  at  the  open  boundaries.  The  solid  line  is  computed  by 
averaging  water  levels  from  every  grid  point  in  the  Gulf  at  each  time  step.  The  dashed  line  is 
computed  from  the  continuity  equation.  The  datum  is  mean  sea  level. 
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simulation  of  HUR23  with  radiation  condition  at  the  open  boundaries.  The  solid  line  is  computed  by 
averaging  water  levels  from  every  grid  point  in  the  Gulf  at  each  time  step.  The  dashed  line  is 
computed  from  the  continuity  equation.  The  datum  is  mean  sea  level. 
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computed  by  averaging  water  levels  from  every  grid  point  in  the  Gulf  at  each  time  step.  The  dashed 
line  is  computed  from  the  continuity  equation.  The  datum  is  mean  sea  level. 


initial  perturbation  of  ? jG  before  the  maximum  peak  surge  is  retained, 
but  is  reduced  on  the  order  of  30  %.  Changes  in  the  maximum  peak 
surge  at  the  coast  are  shown  in  Table  3.  The  magnitudes  of  ??G 
(|j]Gj),  determined  at  the  time  of  the  maximum  peak  at  the  coast,  are 
shown  in  column  5.  The  differences  of  |rjG|  between  each  pair  of 
simulations  are  shown  in  column  6.  The  superscript  plus  or  minus 
indicates  whether  |rjG|  obtained  from  the  simulations  with  the 
radiation  condition  is  larger  or  smaller  with  respect  to  the  former 
result.  The  differences  in  |rjG|  are  clearly  responsible  for  the 
changes  in  maximum  peak  surge  at  the  coast. 

d)  Limited  area  model 

Results  obtained  from  the  simulations  of  the  hurricanes  of 
record  and  the  synthetic  storms  strongly  indicate  close  correlation 
between  the  generation  of  a  simultaneous  Gulf-wide  oscillation  and 
the  net  periodic  volume  transport  through  the  ports.  It  is  of 
interest,  therefore,  to  determine  whether  rjG  still  exists  in  the 
absence  of  wind  stress  in  the  deep  water,  including  Florida  and 
Yucatan  Straits.  Therefore,  simulations  of  HUR5,  HUR23  and  hurricane 
Carla  were  repeated  with  a  modified  wind  field  where  the  wind 
stresses  are  set  to  zero  everywhere  in  the  deep  water  (water  depth 
greater  than  the  upper  layer  thickness).  The  continental  shelf 
regions  of  these  modified  HUR5,  HUR23  and  Carla  resemble  limited  area 
coastal  or  shelf  surge  models.  In  the  folio*  *.ng  discussion  these 
models  are  referred  to  as  limited  area  mo'*  aid  their  designated 
names  are  followed  by  (L). 
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Table  3.  Results  of  the  simulations  of  HUR5,  HUR23  and  hurricane 
Allen  with  and  without  radiation  boundary  condition. 


PEAK 

PEAK 

RATIO  OF 

STORM 

COAST 

nG 

PEAK  nG 

!nGl 

A|ncl 

(m) 

(m) 

(m) 

(m) 

HUR5 

6.50 

0.40 

0.35 

HUR5 (R) 

6.31 

0.28 

0.70 

0. 15 

0 . 20~ 

HUR23 

8.52 

0.23 

0.10 

HUR23 (R) 

8.61 

0.14 

0.74 

0.20 

0.10 

ALLEN 

2.40 

0.35 

-0. 10 

ALLEN (R) 

2.65 

0.22 

0.69 

0.10 

0.20 

The  Gulf-wide  oscillation  is  still  excited  in  the  absence  of 
wind  stresses  in  the  deep  water.  The  time  sequence  of  the  tjg 
response  obtained  from  the  simulation  of  HUR5(L)  is  shown  in  Pig. 

109.  The  maximum  peak  of  this  signal  is  0.2  m  and  the  average  period 
is  28  h.  Magnitudes  of  r?G  obtained  from  the  simulations  of  HUR5  and 
HUR5 ( L )  differ  by  a  factor  of  2.  Neglecting  the  deep-water  wind 
driven  flow  produces  a  dramatic  drop  of  the  total  transport  through 
the  ports  for  HUR5(L)  as  compared  to  HUR5  (not  shown).  The  volume 
transport  through  both  Florida  and  Yucatan  Straits  obtained  from 
HUR5 ( L )  is  caused  by  the  tilting  of  the  water  surface  toward  storm 
center  due  to  the  atmospheric  pressure. 

Simulation  of  HUR23(L)  yields  an  rjQ  response  (Fig.  110)  that  is 
almost  identical  to  that  obtained  from  HUR23.  Comparison  of  the  time 
histories  of  the  volume  transports  obtained  from  the  two  runs  (not 
shown;  reveals  that  the  net  periodic  volume  transports  are  about  the 
same  (6.6  x  106  nrVs)  at  their  first  maxima.  Figure  111  shows  the 
gulf  wide  oscillation  obtained  from  hurricane  Carla(L)  simulation 
which  is  again  comparable  to  that  obtained  from  the  original 
hurricane  Carla  result.  The  initial  positions  of  HUR5  and  hurricane 
Carla  appear  to  be  a  key  factor  in  the  almost  duplicate  nG  responses 
generated  by  the  limited  area  model  and  the  complete  model.  Both 
storms  are  initially  located  closer  and  more  directly  to  the  south  of 
Yucatan  Strait  as  compared  to  HUR5 .  The  stronger  atmospheric 
pressure  gradient  through  Yucatan  channel  is  apparently  more 
effective  in  drawing  water  through  this  opening. 


rricane  5(L) 

M  Hydrograph  (ETAG) 
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lg.  109.  Time  sequences  of  the  average  water  levels  In  the  Gulf  of  Mexico,  r?G,  obtained  from  the 
sicsulati  .  of  HURS(L).  The  solid  line  is  computed  by  averaging  water  levels  from  every  grid  point  i 
the  Gulf  at  each  time  step.  The  dashed  line  is  computed  from  the  continuity  equation.  Th-  datum  is 
mean  sea  level. 
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simulation  cf  hurricane  Carlaib)  .  The  datum  is  mean  sea  level. 


Another  important  result  found  in  the  three  limited  area  model 
simulations  is  the  appreciable  water  surface  elevation  obtained  at 
the  edge  of  the  shelf.  This  water  level  perturbation,  excluding  the 
inverted  barometric  height  at  the  edge  of  the  shelf,  41,  reaches  a 
maximum  of  0.1  m-  0.3  m.  In  the  simulations  of  the  same  storms  with 
the  actual  wind  field,  the  maximum  increases  to  0.7  m  -  0.8  m.  The 
maximum  of  and  the  peak  surge  at  the  coast,  locations  where  the  two 
peaks  occur  and  the  ratio  of  the  two  peaks  obtained  from  all 
simulations  are  presented  in  Table  1. 

It  is  important  to  note  that  the  differences  in  the  maximum  peak 
surges  at  the  coast  between  the  full  and  limited  area  model 
simulations  are  comparable  to  the  differences  in  the  +  maxima.  The 
presence  of  4>  is  important  to  a  properly  posed  open  boundary 
condition  for  limited  area  coastal  surge  models.  Generally,  these 
models  neglect  by  specifying  a  constant  water  level,  which  is 
equivalent  to  the  inverted  barometric  height,  as  the  open  boundary 
condition.  As  noted  above,  this  condition  might  result  in  ar. 
underestimated  maximum  peak  surge  at  the  coast  by  a  value  comparable 
to  the  neglected  ♦.  Another  salient  feature  of  ♦  pertaining  to  the 
consideration  of  open  boundary  conditions  is  the  nonuniforo 
distribution  along  the  shelf  break.  Figs.  112  through  114  show 
profiles  of  ♦  obtained  from  the  original  versions  of  HUR5,  HUR23  and 
hurricane  Carla.  The  nonuniformity  of  these  profiles  implies  that 
the  effect  of  ♦  on  the  peak  surge  at  the  coast  is  different  from  one 
grid  block  to  the  other.  This  result  ultimately  prompted  a  question 
on  specifying  a  constant  water  level  along  the  shelf  break  as  the 
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I  are  the  sanie  as  in  Table  2 

is  the  ratio  between  the  peak  surge  at  the  shelf  break  and  the  peak  surge  at  the  coast. 


ig.  112.  Dynamic  peak  surge,  *,  profile  along  shelf  break  obtained  from  the  simulation  of  HUR5. 
Grid  indicies  are  along  the  bottom  and  the  time  of  peak  *  is  given  by  day  and  hour  across  the  top. 
The  datum  is  mean  sea  level. 
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Fig.  113.  Dynamic  peak  surge,  ♦,  profile  along  shelf  break  obtained  from  the  simulation  of  HUR23 . 
Grid  indicies  are  along  the  bottom  and  the  time  of  peak  ♦  is  given  by  day  and  hour  across  the  top 
The  datum  is  mean  sea  level. 
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ig.  114.  Dynamic  peak  surge,  ♦,  profile  along  shelf  break  obtained  from  the  simulation  of  hurricane 
Carla.  Grid  iodides  are  along  the  bottom  and  the  time  of  peak  *  is  given  by  day  and  hour  across  the 
top.  The  datum  is  mean  sea  level. 
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open  boundary  condition.  The  time  at  which  the  maximum  *  occurs, 
shown  at  the  top  of  the  figures,  introduces  the  time  scale  as  another 
factor  to  be  considered.  In  addition,  the  results  shown  in  Table  4 
indicate  that  the  magnitude  of  '!’  decreases  with  increasing  storm 
forward  speed.  Accordingly,  there  appears  to  be  no  simple  resolution 
to  the  question  of  what  is  the  proper  open  boundary  condition  for 
limited  area  coastal  surge  models.  A  possible  pragmatic  solution  is 
to  use  the  time  sequence  of  water  levels  at  each  grid  point  along  the 
shelf  break  obtained  from  an  a  priori  simulation  covering  the  entire 


basin. 


CHAPTER  V 


SUMMARY  AND  CONCLUSIONS 

A  time-dependent,  numerical,  normal  mode  model  portraying  the 
linear  (except  for  dissipation)  physics  of  a  two  layer  Gulf  of  Mexico 
has  been  developed,  tested  and  verified.  Comparisons  of  the 
numerical  hydrographs  and  the  known  local  responses  produced  by 
hurricanes  Carla  and  Allen  generally  verify  the  barotropic  surge 
response  within  the  Gulf  of  Mexico.  Those  significant  deviations 
which  occur  can  be  explained  by  unique  gauge  locations, 
insufficiently  resolved  model  renditions  of  the  coastline  and 
bathymetry,  and  possible  inaccuracies  in  the  stipulated  atmospheric 
forcing  fields.  The  model  was  used  to  obtain  the  responses  of  the 
Gulf  to  a  series  of  synthetic  storms.  These  responses  provided 
information  which  were  used  to  answer  basic  questions  on  forerunner 
surges  in  the  Gulf  and  related  aspects  of  the  surge  prediction 
problem.  These  questions  concerned  the  influence  of  baroclinic 
motions  on  the  nearshore  surge,  the  establishment  of  forerunners  and 
the  time  and  space  scales  of  this  initial  rise  in  water  level,  and 
surge  conditions  at  the  shelf  break. 

Results  from  the  simulations  of  the  pure  barotropic  and  the  two 
mode  models  of  HUR5  (a  synthetic  storm  with  30  ka  radius  of  maximum 
wind  translating  at  35  km/h  from  the  Cayman  Sea  through  Yucatan 
Strait  and  landfalls  at  Corpus  Christi,  PATHl)  demonstrate  the 
insignificant  contribution  of  the  baroclinic  responses  to  the  water 
levels  on  the  shelf  areas.  The  maximum  surges  in  particular  are 


scarcely  affected  by  the  baroclinic  deep-ocean  responses  in  spite  of 
the  fact  that  such  modes  contain  significant  energy.  In  addition, 
the  initial  rise  of  water  level  before  the  peak  surge  (which,  by 
definition,  is  the  forerunner) in  the  Galveston  hydrographs  obtained 
from  the  two  versions  of  HUR5  show  no  visual  differences.  Based  on 
these  results,  the  baroclinic  response  is  not  important  in  the 
forerunner  surge  phenomenon. 

The  quasi-linear ,  coupled,  normal  mode  model  shows  that  the 
hurricane  induced  forerunner  surge  in  the  Gulf  of  Mexico  is 
associated  with  a  Gulf-wide  oscillation  of  water  level,  tjg.  The 
ubiquity  of  the  rjG  signal  in  the  hydrographs  from  stations  around  the 
Gulf,  except  near  Florida  Strait,  indicates  that  r)G  is  dominated  by  a 
volume  (i.e.,  Helmholtz)  mode.  The  Helmholtz  mode  is  characterized 
by  a  relatively  uniform  amplitude  and  phase,  except  near  the  open 
ports  where  the  amplitude  and  phase  changes  rapidly.  The  forerunner, 
therefore,  has  space  scales  comparable  to  the  horizontal  dimensions 
Of  the  Gulf  of  Mexico. 

Examination  of  the  transport  through  the  ports,  particularly  for 
the  21  day  simulations  for  PATH1,  reveal  large  and  nearly  equal 
amplitude  but  out-of-phase  oscillations  of  about  6.5  day  (156  h) 
period.  Superimposed  on  the  6.5  day  oscillations  are  smaller 
amplitude  in-phase  oscillations  with  a  period  of  about  28  h.  The  out  - 
of -phase  transport  implies  that  when  the  transport  is  in  one  port  it 
is  out  of  the  other.  It  is  the  in-phase  oscillations  in  transport 
which  are  associated  with  the  Helmholtz  mode  in  the  Gulf .  The  long 
period  out-of-phase  oscillations  in  transport,  on  the  other  hand,  are 


associated  with  a  quasi-geostrophic  tilt  of  water  level  across  the 
ports  as  well  as  a  tilt  from  southeast  to  northwest  across  the  Gulf. 
The  6.5  day  period  oscillations  coincide  with  the  natural  modes  on 
the  west  Florida  shelf  (Marmorino,  1982).  The  excitation  of  the  out- 
of-phase  transport  and  associated  quasi-geostrophic  tilt  mode  in  the 
present  model  is  due  to  out-of-phase  wind  forcing  at  the  ports  (as 
can  be  produced  by  the  cyclonic  circulation  in  those  hurricanes  which 
traverse  the  Cayman  Sea) .  Observational  evidence  to  confirm  the 
existence  of  the  6.5  day  mode  in  the  actual  Gulf  is  lacking  (e.g., 
the  Key  West  gauge  records  do  not  show  clear  evidence  for  the  large 
amplitude  6.5  day  oscillations  which  the  model  simulations  for 
hurricane  Allen  or  HUR5  produce).  However,  except  for  locations  near 
the  ports,  the  water  level  variations  associated  with  the  6.5  day 
tilt  mode  are  small  (less  than  0.1  m  along  the  northern  and  western 
coast  of  the  Gulf ) .  Hence  the  quasi-geostrophic  tilt  mode  is 
probably  not  important  with  respect  to  forerunners  in  the  northern 
and  western  parts  of  the  Gulf . 

Close  examination  of  the  ?)q  time  sequences  for  the  long  term 
simulations  for  PATH1  show  that,  in  addition  to  the  presence  of  the 
Helmholtz  mode,  a  period  of  about  3.4  days  is  also  present.  But 
there  is  very  little  evidence  of  the  6.5  day  tilt  mode  in  The 
3.4  day  mode  in  fact  shows  up  in  the  time  history  of  the  spatial  mean 
water  level  tor  the  Cayman  Sea.  Some  observational  evidence  for  such 
a  mode  of  oscillation  exists  for  the  Gulf  (Halper,  1964,  Kelly,  1965, 
Kirvan  et  al.,  1964).  The  relative  excitation  of  these  modes  during 
the  forcing  stage  by  hurricanes  depends  upon  the  path,  the  storm 


scale  and  intensity,  and  to  a  lesser  extent  on  the  translational 
speed  of  the  hurricane.  In  general,  the  tjg  during  the  first  few  days 
after  the  hurricane  enters  the  Gulf  consists  of  three  components. 

The  first  two  components  (i.e.,  the  Helmholtz  mode  and  the  3.4  day 
Gulf-Cayman  mode)  are  associated  with  the  volume  transport  through 
the  ports.  The  third  component  is  a  directly  forced  response 
associated  with  the  spatial  average  value  of  b  (the  inverted 
barometer  term)  over  the  Gulf.  PATH1  tends  to  give  the  largest  peak 
tjg,  particularly  for  larger  radius  storms  (compare  HUR1  and  HUR4 
results) . 

Results  from  the  limited  area  model  and  the  full  model  disclose 
that  both  the  central  pressure  deficit  and  the  wind  induced 
transports  through  the  ports  can  excite  the  Helmholtz  mode.  However, 
the  relative  importance  of  the  two  forcing  fields  in  generating  j}q 
depends  upon  the  storms'  paths  and  their  evolution.  The  almost 
identical  r)G  obtained  from  the  full  and  limited  area  models  of 
hurricane  Carla  and  HUR23  show  that  the  atmospheric  pressure  gradient 
through  Yucatan  Strait  was  more  important  in  the  generation  of  tjg 
than  the  wind  forcing.  In  contrast,  the  small  tjg  response  in  the 
absence  of  wind  in  deep  water  for  HUR5(l)  implies  that  in  this  case 
wind  forcing  was  the  major  factor  in  generating  nG. 

The  average  periods  of  the  Helmholtz  mode  from  this  study  are  in 
the  range  of  25-32  hours.  Platztaan  (1972)  obtained  a  free  Helmholtz 
mode  with  a  period  of  only  21.2  hours.  This  difference  in  periods  is 
due  to  the  fact  that  gG  is  composed  of  both  forced  and  free 
components.  The  average  period  of  26  hours  obtained  from  the  21  day 
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simulation  of  HUR5  was  estimated  during  the  later  stage  of  the 
simulation  for  which  the  forced  component  of  tjg  probably  has  been 
damped  out.  The  variation  of  the  Deriod  of  t?g,  from  one  simulation 
(of  the  original  synthetic  storms  HUR1  through  HUR25)  to  the  other 
may  be  caused  by  the  variation  of  the  period  of  the  forced  component 
which  is  subjected  to  different  forcing.  There  is  some  indication 
that  the  period  of  r)G  is  smaller  for  a  faster  moving  storm.  Based  on 
the  Gulf  mean  tidal  response  in  the  diurnal  band,  the  (JOMT  model  of 
Reid  and  Whitaker  (1981)  used  a  28.5  hour  Helmholtz  mode.  Their 
volume  mode  period  was  adjustable  because  the  radiation  boundary 
condition  employed  in  the  GOMT  model  effectively  takes  into  account 
the  added  mass  of  the  Cayman  Sea  (including  the  Atlantic  Ocean)  by 
means  of  a  complex  ocean  impedance.  The  28.5  hour  period  used  by 
Reid  and  Whitaker  (1981)  is  close  to  the  average  period  of  the  free 
component  of  ijg  obtained  from  this  study. 

It  is  likely  that  Gulf  hurricanes  in  general  elicit  the 
Helmholtz  mode  and  the  longer  period  modes,  but  only  certain  storms 
generate  a  forerunner.  With  tne  definition  of  a  forerunner  as  the 
initial  rise  of  water  level  before  the  peak  surge,  storms  traversing 
PATH4  do  not  generate  forerunner  surges  but  certainly  excite  The 

synthetic  storm  HUR25  (along  PATH5)  serves  as  another  example  of  the 
situation  where  ?}q  exists  but  with  no  forerunner.  The  presence  of  a 
forerunner  in  local  hydrographs  is  therefore  dependent  upon  the  path, 
but  is  also  dependent  on  the  landfall  position  and  the  time  of  the 
excitation  of  nG  relative  to  the  peak  surge.  This  would  explain  why 
every  hurricane  of  record  does  not  have  an  associated  forerunner. 


For  example,  a  fast  moving  hurricane,  which  is  generated  inside  the 
Gulf  might  have  an  associated  initial  peak  of  tjq  which  is  nearly  in 
phase  with  the  primary  shelf  surge. 

A  supplementary  question  addressed  in  the  present  study  concerns 
the  surge  behavior  at  the  shelf  break.  This  is  at  least  as  important 
as  the  forerunner  behavior  for  limited  area  models.  In  many 
applications  of  surge  models  the  domain  of  the  model  is  limited  to  a 
section  of  the  continental  shelf  extending  seaward  from  shore  to  the 
shelf  break  (about  20o  m  depth).  A  common  boundary  condition 
employed  in  such  limited  area  models  is  to  set  the  water  level  (r>) 
equal  to  the  local  value  of  b  for  a  given  time  during  the  traverse  of 
the  hurricane  through  the  model  domain.  The  present  study  in  which 
the  whole  Gulf  of  Mexico  (and  part  of  the  Cayman  Sea)  is  modeled, 
shows  that  the  water  level  at  the  shelf  break  can  depart 
significantly  from  b  at  the  shelf  break.  Moreover,  this  departure 
(n-b)s  has  a  behavior  differing  from  that  of  i)G  and  generally  of 
larger  magnitude.  For  example,  large  scale  *  60  km)  nurricane 

simulations  along  PATH4  yield  larger  values  of  (i)-b)s  than  do  those 
along  PATH1 ;  this  is  just  the  opposite  behavior  of  tjq  for  these  two 
paths. 

In  order  to  gain  some  further  insight  with  respect  to  the  shelf 
break  condition,  three  hurricane  runs  (HUR5,  HUR23  and  Carla)  were 
repeated  with  the  winds  turned  off  in  the  deep  region  of  the  Gulf  and 
Cayman  Sea  (i.e.,  for  depths  greater  than  H^) .  This  is  equivalent  to 
having  a  limited  area  shelf  model  (including  all  shelves  in  the 
system)  but  allowing  wave  energy  to  radiate  into  the  deep  Gulf. 
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Comparison  of  the  tjs  (i.e.,at  the  shelf  break)  from  these  runs  with 
their  counterparts  for  the  fully  forced  model  shows  that  there 
remains  a  significant  difference  (of  order  of  0.4  m)  in  the  peak 
values. 

The  primary  conclusions  of  this  study  can  be  summarized  as 
follows: 

(1)  The  surge  on  the  shelf  including  the  forerunner 

is  primarily  a  barotrpic  response;  very  little  of  the 
berociinic  energy  generated  within  the  deep  water  regions  is 
transmitted  onto  the  shelf. 

(2)  The  forerunner,  when  it  occurs,  is  associated  primarily 
with  the  Gulf -wide  modes  contained  in  ?}q  (the  spatial 
average  of  the  Gulf  water  level  at  a  given  time);  this  is  a 
volume  mode  which  exhibits  periods  of  about  28  hours  and  3.4 
days. 

(3)  Forerunners  are  always  associated  with  »)G,  but  not  all 
hurricanes  which  excite  have  an  associated  forerunner. 
Regardless  of  whether  or  not  forerunners  exist,  jjq  can,  if 
properly  phased,  affect  the  magnitude  of  the  peak  surge. 

(4)  Limited  area  shelf  models  (at  least  within  the  Gulf  of 
Mexico),  which  employ  the  seaward  boundary  condition  ty  ■  b 
(or  the  generalization  of  this  which  allows  outward 
radiation  of  free  waves)  will  always  underestimate  the  peak 
surge  at  shore;  the  underestimate  can  amount  to  as  much  as 


10  percent. 


(5)  A  Gulf-wida  quasi-geostrophic  tilt  mode  of  about  6.5  day 
period  is  found  in  the  model  superimposed  or  the  volume 
mode,  but  verification  of  this  from  observations  is  lacking. 

The  above  conclusions  are  based  on  a  quasi-linear,  two-layer 
model  of  the  Gulf  of  Mexico  and  a  portion  of  the  Cayman  Sea.  The 
model  used  a  grid  size  of  15'  in  latitude  and  longitude  and  allowed 
for  a  variable  Coriolis  parameter.  The  lack  of  non-linear  advection 
and  the  ad  hoc  conditions  at  the  open  boundaries  of  the  Cayman  Sea 
and  Florida  Strait  should  be  borne  in  mind,  particularly  with  regard 
to  conclusion  (5).  Long  period,  quasi-geastrophic  (planetary)  modes 
are  known  to  exist  in  the  ocean;  however,  their  spatial  structure 
and  behavior  is  known  to  be  sensitive  to  open  boundary  conditions  and 
to  non-linear  phenomena  such  as  advection  of  vorticity  (which  is  not 
admitted  in  the  present  model).  Indeed,  the  effect  of  the  strong 
quasi-steady  Loop  Current  within  the  eastern  Gulf  is  missing  in  the 
present  linear  model.  Regardless  of  these  limitations,  it  is  felt 
that  conclusions  (1)  through  (4)  remain  valid. 
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